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INFLUÊNCIA DAS TENSÕES NEUTRAS NEGATIVAS 
NAS CARACTERÍSTICAS ESTRUTURAIS 
DOS SOLOS COMPACTADOS 


Conclusão) 


As curvas representadas nas figuras 4.3 e 4,4 


traduzem as variações da componente distorcio- 


nal do estado de tensão efectiva (0 ), da tensão 


Clbhgiem) 
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Fig. 4.3 — Compressão axial de amostras 
não saturadas (1.º série) 
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neutra na fase líquida (u, —u,,) e da deformação 


volumétrica relativa (:,), em função da extensão 


axial (5,). 

O comportamento geral na 1.º série caracteri- 
za-se inicialmente por uma diminuição em valor 
absoluto das pressões neutras na fase líquida, 
facto que se conjuga bem com a diminuição do 
indice de vazios registada. A resposta negativa- 
mente dilatante nesta fase inicial é a esperada 
dum ensaio drenado, dado o efeito da elevada 
taxa de aumento de q, na tensão média. Mas, à 
medida que se processou a deformação axial, 
verificou-se uma inversão no sentido da variação 
das deformações volumétricas. Na realidade a 
dimiminuição inicial e em valor absoluto das 
pressões neutras negativas repercutiu-se na 
pressão média efectiva, fenómeno que, associado 
a uma tensão distorcional crescente, determinou 
aumento do índice de vazios. 

Quanto aos resultados dos ensaios da segunda 
serie (figura 4.4), verifica-se que, quando compa- 
rados com os da primeira série, apresentam as 
seguintes diferenças: 

a) Os valores máximos de 3, registaram-se para 
maiores deformações axiais. 

b) As tensões neutras sofreram variações de 
maior amplitude. Mantendo-se no campo nega- 
tivo, a evolução deu-se sempre no sentido duma 
diminuição em valor absoluto e com um anda- 
mento do mesmo tipo para todas as amostras. 

c) As inversões de sentido das variações volu- 
métricas registaram-se para deformações axiais 
tanto maiores quanto maior o grau de saturação 
inicial. 

d) Os valores máximos de c1, em função do 
grau de saturação inicial, obtidos para as duas 
séries, estão representados na figura 4.5. Com- 
parando as duas curvas verifica-se que, para o 
mesmo grau de saturação inicial, a resistência 
ao corte é menor quando o ensaio se processou 
a velocidade de deformação maior. 
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Hã assim indicações de que os ensaios da pri- 
meira série foram demasiado rápidos. Na reali- 
dade, como o tempo de repouso dos provetes 
não afectou de maneira sensivel as pressões neu- 
tras negativas iniciais (figura 4.2) é de crer que 
para a primeira série as tensões neutras na zona 
de corte, embora negativas, tenham atingido 
menores valores absolutos que os registados. Os 
ensaios não foram suficientemente lentos para 
evitar diferenças acentuadas entre as pressões 
neutras instaladas ma zona média da amostra e 
as lidas na base. 

Analisados os resultados dos ensaios de 
compressão axial de amostras não saturadas 
conclui-se que: 


1) A máxima resistência ao corte depende da 
velocidade de deformação, sendo, para velo- 
cidades mais lentas, menor essa resistência, 


2) Nos ensaios com menor velocidade de 
deformação, a máxima resistência ao corte 
estã associada a uma maior deformação. 


3) A medição das tensões neutras é também 
dependente da velocidade de deformação. 
Como a relação tensão-deformação é influen- 
ciada pela tensão instalada na fase líquida, 
pode julgar-se da adequabilidade da velo- 
cidade adoptada analisando as variações 
volumétricas da amostra. Assim, se o indice 
de vazios mínimo coincidia com , máximo, 
a pressão neutra registada deve traduzir a 
pressão neutra equivalente realmente insta- 
lada. Os ensaios da segunda série obedecem 
sensivelmente a esta condição. 


4) Como é bem sabido, o comportamento, 
quanto a deformações volumétricas, dum 
solo normalmente consolidado e dum solo 
pré-consolidado em ensaio drenado, é o 
representado na figura 4.7. Analisando as 
deformações volumétricas dos ensaios reali- 
zados (figuras 4.3 e 4,4) verifica-se estar 
em presença de solos pré-consolidados. 
Como seria de esperar, o grau de pré- 
-consolidação é uma função do grau de 
saturação inicial. Pretendeu-se pois verificar 
se esse efeito, presente no solo não saturado, 
se mantinha após a saturação, condicionando 
assim o comportamento do solo em termos 
de resistência ao corte. 
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4,2 — Efeitos de pré-consolidação, em solos 
compactados saturados, resultantes de 
tensões neutras negativas previamente 
instaladas na fase liquida 


Com as amostras 32, 33 e 35 realizaram-se 
ensaios de corte em compressão triaxial do tipo 
consolidado não drenado. Para assegurar a satu- 
ração utilizou-se uma contrapressão de 8,0 kg/cmº. 
Os resultados obtidos apresentam-se na figura 4.6 
sendo 7º = 33º. Admitindo que, para um dado 
indice de vazios inicial, os valores de coesão (c') 
e ângulo de atrito (1') são independentes do grau 
de saturação (Bishop, 1960) comparam-se, na 
situação de rotura, os resultados dos ensaios com 
solos saturados e não saturados da 1.º série 
(figura 4.8), sendo que nestes últimos se admitiu 
que a rotura se processou aquando da inversão no 
sentido de variação das deformações volumétricas. 


Fig. 4.4 — Compressão axial de amostras 
(não saturadas 2.* série) 
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Obtiveram-se as curvas 1 e 2 respectivamente 
para as amostras não saturadas e saturadas. De 
acordo com os círculos de Mchr representados 
na figura, o designado por | traduz uma condição 
em que a razão ZM, (us, —u,), tem um valor 
inferior à unidade. Mas, dado que a curva 1 
intersecta o eixo das ordenadas, para uma amostra 
correspondente a circulo II o valor obtido para 


et — = i F n 1 ” 
a relação 24! (u, — Vw Im é superior à unidade. 


Isto significa que a saturação, ao anular as tensões 
neutras negativas (portanto a correspondente 
tensão efectiva), não implica o total desapare- 
cimento duma tensão efectiva residual. Por outras 
palavras, a curva 2 não deve passar na origem 
mas sim intersectar, com valor positivo, o eixo 
das ordenadas. Pode explicar-se este comporta- 
mento como sendo o efeito de pré-consolidação 
relacionado com as tensões neutras negativas 
previamente instaladas e dependente, em primeira 
aproximação, da ordem de grandeza e tempo de 
actuação dessas tensões. 

Do ponto de vista prático pretende-se escla- 
recer a questão, ainda hoje controversa, de se 
considerar ou não uma «parcela coesiva» na 
resistência ao corte mobilizável, após saturação, 
em aterros argilosos compactados. 

Com amostras da segunda série realizaram-se 
ensaios triaxiais em que, de acordo com os resul- 
tados dos ensaios anteriores, se utilizaram pressões 
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Fig. 4.5 — Amostras não saturadas, Variação 
de «| com o grau de saturação de fabrico 
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de confinamento substancialmente inferiores 
(Quadro 4.II). Também a saturação envolveu 
cuidados especiais, pois pretendendo-se graus de 
saturação entre 90 a 100º/o ter-se-ia de recorrer 
a contrapressões elevadas. Foi o caso dos ensaios 
com as amostras da primeira série. 

Mas se o emprego de contrapressão dessa 
ordem de grandeza não levanta problemas nos 
solos arenosos, podendo considerar-se, pelo menos 
dum ponto de vista prático, como não interferindo 
no estado de tensão efectiva instalado, outro 
tanto pode não suceder com os solos argilosos. É 
uma questão estreitamente vinculada à proble- 
mática que relaciona áreas de contacto com 
tensões efectivas, questão ainda carecida de 
investigação. 

Como ao saturar se pretendeu fundamental- 
mente eliminar a interferência da tensão superfi- 
cial mas tensões entre partículas, justifica-se em 
seguida como se pode obter o mesmo resultado 
sem forçar a saturação total. 
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Fig. 4.6 — Ensaios triaxiais (consolidados, não drenados 
e saturados com uma contrapressão de 8,0 kg/em?) 


Na fig. 4.9 representa-se em esquema um 
plano ondulado contendo as áreas de contacto 
entre as partículas dum solo não saturado. 

O estado de tensão dum elemento de solo 
pode definir-se a partir do tensor q, e das pre- 
sões neutras u, e u,. O tensor das tensões 
totais, 7, resulta da sobreposição dum tensor 
esférico, p, e dum distorsional 9. Assim, 

=P dj + (4.1) 
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Fig. 4.7 — Deformações volumétricas 
nos ensaios 
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em que 9, é o simbolo de Kronecker. De acordo 
com o esquema da figura 4.9 pode estabelecer-se 
a equação de equilíbrio 


di Tae e Ay Uw “ij As e 


l 
—dy JT dx (4.2) 
em que 
T ==tensão superficial por unidade de compri- 
mento dos meniscos 
x == comprimento da intersecção da interface 
ar-água com o plano 
ci =tensão equivalente à transmitida pelos 
contactos sólidos 


A = Zica: A, = Dom; a=a,taçta, 


w 
f T dx == termo que diz respeito à força (por área 
unitária do plano) que é induzida pela 
tensão superficial da interface ar-água 
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Fig. 4.8 — Ensaios de amostras saturadas 
e não saturadas (1.º série) 


Combinando as equações 4.1 e 4.2 e admitindo 
que A, é desprezável relativamente a À e A, 
(o que só não será aceitável para argilas com 
disposição lamelar das partículas), obtém-se 


C,; =|p—u, + A, (u,—u,) + [T dx [>< 


X 9 + , (4.3) 
Quando da aplicação de uma tensão hidrostá- 
tica a tensão intergranular será 


c=ec—-u+A, (us —u,) + | T dx (4.4) 


sendo portanto função duma componente exte- 
rior (7 —u,) e de duas internas—(u,-u,) e 
ST dx. Existe pois uma diferença importante 
entre estas duas últimas componentes do estado 
de tensão. Enquanto (u, —u,) actua uniforme- 
mente sobre toda a superfície molhada das par- 
tículas, as forças devidas à tensão superficial 
concentram-se nas pequenas áreas dos contactos 
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sólidos. Então, uma variação na interface ar-água 
não só produziria uma alteração na tensão inter- 
granular como poderia originar uma redistribui- 
ção na tensão entre partículas (Matyas e Radha- 
krishna, 1968). Assim, nas amostras ensaiadas 
procurou-se saturar até um grau a que corres- 
pondesse uma presença de ar sob a forma de 
bolhas esféricas o que, do ponto de vista em 
questão, é análogo à saturação total. Na reali- 
dade a tensão | T dx anula-se e da tensão 
(u, —u,) mediu-se a parcela correspondente a 
u,- À pressão u, no interior das bolhas não é 
conhecida sendo que, para o presente caso, tal 
facto não apresenta inconvenientes. 

As amostras foram inicialmente consolidadas 
a teor de humidade constante. Seguidamente foi 
facultado o acesso de água ao provete e, uma vez 
estabelecido equilíbrio, iniciou-se a aplicação lenta 
de contrapressão até atingir um grau de satu- 
ração de cerca de 96%/,. A técnica deste ensaio 
já foi referida em 3, devendo destacar-se como 
principais objectivos os seguintes : 


a) Conhecimento do grau de saturação em 
qualquer instante. 

b) História de tensões análoga para os provetes 
dos vários grupos submetidos à mesma tensão 
de consolidação. 


Fig. 4.9 — Modelo de um solo parcialmente saturado. 


Verificou-se que as contrapressões aplicadas 
não ultrapassaram 1,6 kg/cm*”, não estando por- 
tanto as amostras sujeitas a contrapressões 
muito elevadas, que seriam necessárias para obter 
a saturação completa e que, em última análise, 
teriam influência análoga à resultante do em- 
prego de grandes pressões de consolidação. 

A via mais directa para verificar os efeitos da 
pré-consolidação consistiria na análise das envol- 
ventes, na rotura, para cada grupo de amostras. 
Como a diferença entre os grupos (Quadro 4.II) 
é o teor de humidade de fabrico, procura-se 
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assim relacionar esses graus de saturação iniciais 
com resistências ao corte exibidas pelas amos- 
tras após saturação. 

No entanto, como é bem sabido, nos ensaios 
triaxiais consolidados não drenados a tensão de 
confinamento constante os solos normalmente 
consolidados tendem a passar a pré-consolidados. 
Este efeito será tanto mais importante quanto 
maior o aumento da tensão neutra vinculado 
à aplicação do tensor distorcional. 

Para contornar a dificuldade originada pela 
sobreposição nos resultados dos efeitos da pré- 
-consolidação devida às tensões neutras negativas 
e da pré-consolidação possivelmente induzida 
durante o ensaio, reçorreu-se às relações mate- 
máticas que existem entre a envolvente de resis- 
tência ao corte e a função 

c,/7, = F(R,) (4,5) 
(Barron, 1960) em que R, (Hirschfeld, 1958) 
é igual ao quociente entre a tensão de consoli- 
dação 7, e a tensão efectiva normal ao plano de 


rotura c, e se designa por razão de pré-conso- 
lidação. 
Assim paraR, =10 valor de f (R,) éf(1) e 


EO =... 
(MD) +1 
em que vw representa a inclinação da envolvente 


relativa a ensaios consolidados drenados condu- 
zidos sobre amostras do mesmo solo, Mas 


sen (4.6) 


T. ç, + CT; | — 0; 
— = ————— —  — 6 é Pá 4,7 
R, 7 + tor (47) 


DS 


QUADRO a. 


Amostras saturadas. Ensaios tipo consolidado-não drenado 
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QUADRO 4.IV 


Amostras saturadas Ensaios tipo consolidado-drenado 


PESO ESPECIF| TEOR DE [GRAU DE SAT|TENSÃO DE |VELOCIDADE DE 
APARENTE SECO| HUMIDADE |NO INÍCIO 00 CONSOLIDAÇÃO DEFONMAÇÃO 


19,26 
18,88 


Combinando as equações 4.5 e 4.6 tem-se que 


E ar AE (4.8) 
' * | HR) e EMO) 
e a tensão de corte, 7, no plano de corte 
— [F(R,)— E cuca 
T= O; ER 71 vVE(I) (4,9) 


Substituindo o valor de 0, dado pela equação 
4.8 vem 


7 ER) es) aa 
— F(Rç)—1 VE) 
C, f(1) + 1 


(4.10) 


obtendo-se assim o valor de 0 ns no plano de ro- 
tura, em função de R,. 

No Quadro 4.IIIl registaram-se os resultados 
dos ensaios das amostras dos três grupos e, como 
havia necessidade de conhecer o valor de 2”, rea- 
lizaram-se ensaios triaxiais do tipo consolidado- 
-drenado com três amostras saturadas. No Qua- 
dro 4.IV e figura 4.10 representam-se os elemen- 
tos de maior interesse relativamente àquelas 
amostras e respectivos ensaios. Com o valor de 7” 
assim obtido (28,5º) calculou-se o valor de f(R ,) 
para R, = 1.0 em 2,82. 

Na figura 4.11 estão traçadas as curvas ima- 
gem da função 7,/7, = f(R,) para cada um dos 
grupos de amostras. Essas mesmas curvas foram 
transferidas para a figura 4.12 através da equa- 
ção 4.10 e para uma pressão de consolidação 
igual à unidade. Como se pode verificar, as cur- 
vas situam-se acima da envolvente de resistência 
ao corte relativa a provetes normalmente conso- 
lidados e estão representadas em termos de ten- 
sões relativas pelo que, admitindo outro valor 
para a tensão de consolidação, resultariam curvas 


que se sobrepunham às já existentes. Tal facto é 
evidente se se atender às expressões deduzidas, 
embora na prática essa coincidência possa não 
ser verificada, dada a diminuição de índice de 
vazios para valores crescentes da tensão de 
consolidação. 
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Fig. 4.10 — Amostras saturadas. Ensaios do tipo 
consolidado-drenado 


Como já se referiu, todas as amostras foram 
fabricadas de molde a exibirem o mesmo peso 
específico e o mesmo tipo de estrutura e, durante 
os ensaios, procurou-se seguir idênticas trajec- 
tórias de tensões. A diferença entre os três grupos 
de amostras diz respeito ao teor de humidade de 
compactação (Quadro 4.11) sendo pois possível 
relacionar cada uma das envolventes obtidas com 
uma dada pressão neutra negativa inicial. 


Gi 


« 


a 


Hig. 4.11 — Relação entre s,/0;e R, para os 
três grupos de amostras ensaiadas 
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Na figura 4,12 a parte das curvas representada 
a traço interrompido significa que não tem apoio 
experimental para valores de 7,/7. inferiores 
a 0.5. No entanto, conhecido o andamento com 
base nos ensaios, parece legítimo admitir que os 
materiais têm uma resposta bem representada, 
do ponto de vista prático, pelas rectas a traço- 
-ponto da figura 4.12. Qualquer destes segmentos, 
mantendo um ângulo 7 com o eixo das abcissas, 
divide aproximadamente em partes iguais a área 
compreendida entre a envolvente relativa ao 
grupo de amostras a que essa recta se refere. 


7 ea 
te U,s -— 


Fig. 4.12 — Envolventes de resistência em termos 
de tensões relativas 


Trata-se portanto duma interpretação em termos 
de valores médios, ditada por razões de ordem 
prática e só aplicável quando as tensões normais 
efectivas no plano de rotura forem inferiores às 
tensões de pré-consolidação. 

Com os resultados obtidos do estudo do efeito 
da pré-consolidação relacionou-se o grau de 
saturação inicial do solo compactado e coesão 
exibida após saturação. Na figura 4.13 apresenta-se 
a curva que traduz essa interdependência. As 
coesões estão expressas em percentagem da tensão 
de pré-consolidação e a análise da curva obtida 
indica que a coesão tende para um valor limite 
quando o grau de saturação inicial diminui. Esta 
verificação confirma os resultados de estudos 
experimentais até hoje efectuados e que visam o 
conhecimento da relação entre a tensão neutra 
negativa e a tensão efectiva dela resultante. 

Quando o grau de saturação inicial ultrapassa 
85"/y verifica-se, para o solo em questão, que a 
coesão residual é nula. Na realidade a tais graus 
de saturação corresponde a formação de vazios 
de ar com configuração de bolhas esféricas. Assim, 


. TECNICA Nº 407 


os efeitos das tensões neutras da fase líquida 
sobre a estrutura sólida são praticamente nulos, 
pelo que o solo não exibe coesão residual, 


Ca) 
Sr, 19/40) 
Fig. 4.13 — Relação entre o grau de saturação 
inicial e as coesões em termos de tensões 
relativas) exibidas pelo solo saturado 


5—- INTERPRETAÇÃO CONJUNTA DOS 
RESULTADOS 


O estudo experimental efectuado abrange o 
solo quer no estado não saturado quer após a 
evolução para um estado saturado. Pretendeu-se 
caracterizar bem o material na situação inicial 
ou de fabrico, para poder determinar a influência 
no comportamento do mesmo material no estado 
final, ou seja, após saturação. As amostras foram 
preparadas de molde a exibirem o mesmo peso 
específico aparente seco, variando apenas o grau 
de saturação. Compactando estaticamente e para 
teores de humidade correspondentes ao lado seco, 
procurou-se minimizar a dispersão das caracte- 
rísticas estruturais. 

Para o solo ensaiado verificou-se que as tensões 
neutras negativas instaladas na fase líquida das 
amostras compactadas dependem fundamental- 
mente do grau de saturação. 

A determinação das tensões neutras negativas 
iniciais foi realizada em duas séries de amostras 
submetidas a compactação estática. Na primeira 
efectuaram-se medições uma semana após a 
compactação ; na segunda, o intervalo foi de duas 
semanas. Embora se tenha anotado uma pequena 
diferença nos resultados (figura 4.2), a dispersão 
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com que estes são medidos não permite atribuir 
grande significado à informação colhida. 

As duas série de amostras não saturadas foram 
submetidas a ensaios de compressão axial com 
medição de tensões neutras e de deformação 
volumétrica. 

Verificou-se que quanto menor o grau de satu- 
ração maior a resistência dos provetes e menor 
a sua deformabilidade. Quanto à velocidade de 
deformação imposta, dado que a segunda série 
foi ensaiada a uma velocidade menor, verificou-se 
que, para o mesmo grau de saturação, as amos- 
tras submetidas a um ensaio mais lento exibiram 
maiores deformações e resistências. Deve no 
entanto salientar-se que a influência da veloci- 
dade de deformação vai sendo atenuada com a 
diminuição do grau de saturação. 

Estes resultados experimentais evidenciam a 
importância da contribuição da tensão superficial 
no comportamento do solo não saturado. Como 
durante o ensaio os vazios de ar intercomunicam, 
a variação de pressão no ar não determina 
variações da tensão efectiva. 


Esta, como já foi referido, é dada pela 
expressão 
=| yu, + As (uam Uy) + | T dx | >< 
x 3, + G". 


onde os símbolos têm significado conhecido. 
Assim, a contribuição de carácter interno para as 
tensões efectivas é 


À 


w 


u,—(1—-AÇ) u, + | Tex 


em que u, e u, não determinam tensões inter- 
particulares em função da sua própria grandeza 
mas sim da área À, e da diferença (u, — u,). 
A contribuição da parcela / Tdx caracteri- 
za-se ainda por um efeito estabilizador do ponto 
de vista da manutenção duma dada estrutura. 
Como no início do ensaio e à medida que se 
vão instalando as tensões tangenciais, o grau de 
saturação aumenta, o valor de / Tdx será fun- 
ção da dimensão do perímetro dos meniscos 
ar-água. Quer isto dizer que, partindo de graus 
de saturação baixos, a contribuição da tensão 
superficial para as tensões entre as partículas 
aumentará com o aumento do perímetro atrás 
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referido. Por sua vez a curvatura da interface 
ar-água diminui, o que implica a diminuição de 
u,—u,. Mas como À , aumenta, as implicações 
desta alteração da tensão efectiva tendem a ser 
atenuadas. 

No conjunto haverá assim uma diminuição na 
contribuição das tensões de origem interna para 
as tensões efectivas. Mas dado o efeito crescente 
e estabilizador de tensão superficial, as tensões 
distorcionais solicitantes depararão com apreciá- 
vel acréscimo de resistências estruturais. No 
entanto deve ter-se presente que o esquema apre- 
sentado não é tão simples como se pode depreen- 
der destas considerações. Assim, devido às soli- 
citações tangenciais, a estrutura que inicialmente 
será do tipo floculado tenderá a passar progres- 
sivamente ao tipo disperso. Terá então de consi- 
derar-se a sobreposição de dois fenómenos, um 
dizendo respeito a variação do índice de vazios e 
outro à alteração da configuração geométrica 
desses mesmos vazios. 

Conclui-se portanto que o comportamento dos 
provetes não saturados, em termos de tensão- 
-deformação, é determinado pelo grau de satura- 
ção de fabrico. A medida que esse valor vai 
sendo menor, maior será a tendência do corpo 
de prova para exibir um andamento do tipo 
«rigido», isto é, nota-se o aparecimento dum 
valor crescente de «peak» da tensão, para valo- 
res sucessivamente decrescentes da deformação. 
Todas as referências feitas à evolução e efeitos 
da tensão superficial se harmonizam com o 
comportamento observado. 

Resta analisar a influência na relação tensão- 
-deformação da velocidade de deformação im- 
posta. Como se verificou, da diminuição do grau 
de saturação inicial resultou que essa influência 
era cada vez menor. De acordo com o anterior- 
mente exposto, a tensão superficial em solos 
pouco saturados desempenharia um papel esta- 
bilizador da estrutura, o qual se acentuaria 
durante a deformação. Para graus de saturação 
mais elevedos, a contribuição da tensão superfi- 
cial seria já independente da deformação, sendo 
pois as resistências tangenciais mobilizadas fun- 
ção do grau de alteração estrutural, Esse grau de 
alteração seria incrementado com o aumento da 
velocidade de deformação. Assim, a possibilidade 
de a tensão superficial se opor a este efeito iria 
diminuindo progressivamente com a elevação do 
grau de saturação de fabrico. 
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Analisado o comportamento do solo não satu- 
rado, de acordo com os resultados dos ensaios a 
que foi submetido, pretendeu-se seguidamente 
determinar a influência que o grau de saturação 
de fabrico poderia ter ma resistência ao corte 
exibida pelo mesmo solo após saturação. Os 
ensaios de corte em compressão triaxial reali- 
zados evidenciaram um efeito de pré-consolidação 
tanto mais acentuado quanto menor o grau de 
saturação inicial dos provetes. Tem-se assim que 
todo o sistema de tensões internas instaladas na 
fase sólida dum solo não saturado condiciona, 
uma vez este saturado e para uma dada gama de 
tensões aplicadas, a resposta a solicitações tan- 
genciais. 

De acordo com as determinações experimentais 
efectuadas, as pressõesneutras negativas instaladas 
na fase líquida do solo compactado atingem valores 
elevados e produzem alterações de ordem estru- 
tural. Aliás, parece importante sublinhar que por 
compactação estática a vários teores de humidade 
seria impossível obter provetes com estruturas 
idênticas. A intenção foi sobretudo, e como já 
foi acentuado, minimizar a dispersão dessas 
características estruturais. Assim, as amostras não 
têm o seu comportamento totalmente determi- 
nado pelo efeito das tensões neutras negativas. 
Tem de atribuir-se uma certa influência ao pró- 
prio fabrico, o que, dum ponto de vista prático, 
não tem qualquer repercussão. Outro aspecto que 
convém assinalar refere-se à trajectória de tensões 
quando se procede à saturação. De acordo com 
os conhecimentos actuais, quando um solo com 
baixo grau de saturação e tendo instalada uma 
dada solicitação exterior é levado à saturação 
pode sofrer um colapso estrutural. Este fenómeno 
é estreitamente dependente do valor da solici- 
tação aplicada. Na realidade, à medida que se 
processa a saturação há uma diminuição de ten- 
são efectiva de carácter interno. Se a solicitação 
exterior aplicada for superior a um certo valor, 
pode dar-se uma instabilização estrutural assina- 
lada por uma rápida e elevada variação volumé- 
trica. 

No caso presente dado que se elevou o grau 
de saturação das amostras para cerca de 96”/, 
(sob tensões hidrostáticas aplicadas de valores 
compreendidos entre 0,165 e 0,660 kg/cm”), houve 
que verificar a ocorrência de possíveis colapsos. 
Na realidade, para o grau de saturação de 96º/,, a 
tensão superficial deixa de ter qualquer contri- 
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buição para a tensão efectiva ou qualquer função 
de estabilização na estrutura. No entanto não foi 
observado qualquer colapso, sendo o andamento 
das variações volumétricas, durante o processo 
de saturação, do tipo do apresentado como 
exemplo para a amostra 55 e que se representa 
na figura 5.1. Portanto as tensões de consolida- 
ção não foram suficientemente elevadas para 
causar qualquer alteração estrutural brusca. 
Finalmente, determinou-se a relação entre o 
grau de saturação inicial (ou de compactação) e 
as coesões exibidas pelo solo após saturação. 
Para o solo estudado verifica-se que as coesões 
variam de maneira sensível para graus de satu- 
ração compreendidos entre 50 e 85º/o. Para graus 
de saturação inicial superiores a 85º/o, a coesão 
no solo saturado é nula. Quando o grau de 
saturação é inferior a 50º“/, a coesão tende para 
um limite máximo que, no caso presente, é de 
cerca de 15º/, da tensão de pré-consolidação. 


1,000 


AMOSTRA 55 


0,800 


CONTRA-PRESSÃO (kg/cm?) 


0,600 


0,400 


0,200 


PRESSÃO DE CONSOLIDAÇÃO « 0,450 (Ng/em?) 


Avelem?) 


-20 


Fig. 5.1 — Consolidação e saturação duma 
amostra para ensaio triaxial 


6 — CONCLUSÕES 


Do estudo efectuado e para o tipo de solo 
ensaiado pode concluir-se que: 


a) As tensões neutras negativas instaladas na 
fase líquida do solo compactado dependem fun- 
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damentalmente do grau de saturação. O equipa- 
mento construído para medição dessas pressões 
revelou-se satisfatório. De acordo com os teores 
em água empregados na preparação das amos- 
tras, mediram-se tensões neutras com valores 
compreendidos entre -1,0 e -8,4 kg/cm*. 


b) Determinaram-se tensões neutras iniciais em 
duas séries de amostras fabricadas por compac- 
tação estática. Na primeira, as medições efec- 
tuaram-se uma semana após a compactação e na 
segunda duas semanas depois. Durante o período 
de repouso as amostras permaneceram conve- 
nientemente isoladas. Verificou-se que a relação 
entre o grau de saturação de fabrico e a tensão 
neutra negativa inicial é idêntica para as duas 
séries; portanto pode inferir-se que, após o pe- 
ríodo inicial em que se processam redistribuições 
de humidade, intervalos de tempo da ordem dos 
que foram ensaiados não têm repercussão signi- 
ficativa nas tensões neutras instaladas. 


c) Amostras parcialmente saturadas ensaiadas 
em compressão axial com medição de tensões 
neutras e de deformações volumétricas eviden- 
ciaram que quanto menor é o grau de saturação 
de fabrico maior é a resistência e menor a defor- 
mabilidade. 

Quanto à influência da velocidade de defor- 
mação imposta, verificou-se que a valores cres- 
centes desta grandeza correspondem valores de- 
crescentes da resistência e deformabilidade. No 
entanto a influência da velocidade de deforma- 
ção tende a ser atenuada com a diminuição do 
grau de saturação inicial. 

Este comportamento deve-se aos efeitos da ten- 
são superficial que actua como elemento condi- 
cionador dominante da estrutura de fabrico e da 
alteração dessa mesma estrutura quando o solo 
sofre deformações relacionadas com uma dada 
solicitação. 


d) O grau de saturação na compactação do solo 
condiciona o seu comportamento quando, após 
saturação, ele é submetido a solicitações tangen- 
ciais. Os resultados dos ensaios de corte com 
compressão triaxial mostraram um efeito de pré- 
-consolidação tanto mais acentuado quanto menor 
o grau de saturação inicial. 

Este comportamento resulta dos efeitos do sis- 
tema interno de forças interpartiçulares, instalado 
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no solo parcialmente saturado e derivado da ten- 
são superficial das interfaces ar-água. 


e) Para o solo estudado quantificou-se o efeito 
do grau de saturação de compactação, determi- 
nando a coesão exibida pelo material após a 
saturação. 

Sendo um facto comprovado que 6 grau de sa- 
turação condiciona a estrutura do solo quando 
da compactação, pode acrescentar-se que esse 
condicionamento estrutural só tem repercussão 
nas resistências ao corte exibidas pelo solo satu- 
rado desde que seja inferior a um dado valor 
crítico. Para valores sucessivamente decrescentes 
em relação a esse valor crítico a influência do 
grau de saturação de fabrico é crescente. Tende 
no entanto para um limite máximo, ao qual cor- 
responde um segundo valor crítico do grau de 
saturação de compactação. Quer isto dizer que 
a estrutura atingiu uma fase singular do seu 
arranjo que corresponde, para as condições de 
compactação, ao máximo grau de floculação. 

Assim, quando se provoca a saturação, quer 
do solo compactado com um grau de saturação 
igual ao valor crítico inferior, quer do solo com- 
pactado com um grau de saturação menor do que 
aquele, como o estado estrutural é idêntico, as 
trajectórias de tensões derivadas do processo 
de saturacão têm repercussões estruturais seme- 
lhantes. Determinam assim a mesma pré-conso- 
lidação no estado final, 
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ANÁLISE DE ESTRUTURAS REDUTIVEIS 
A OSCILADORES NÃO LINEARES COM HISTERESE ! 


RESUMO 


Nos estudos sobre Engenharia Sismica que vêm sendo 
conduzidos no LNEC desde há cerca de 15 anos, tem sido 
necessário aplicar técnicas de análise numérica a osciladores 
de um e vários graus de liberdade, lineares e não lineares, 
os quais constituem os modelos matemátices de diversos tipos 
de estruturas, 

Descrevem-se neste trabalho as características dos osci- 
tadores que têm sido analisados, relacionando-as com os 
tipos de estruturas que representam, Referem-se os princi- 
pais aspectos da formulação das equações diferenciais que 
descrevem o comportamento dos sistemas e das técnicas uti - 
lizadas na respectiva integração «passo a passo», Apre- 
sentam-se vplicações ao estudo do comportamento de edifi- 
cios de betão armado actuados por sismos de forte intensi- 
dade, 


1 — INTRODUÇÃO 


Sob a acção dos sismos, as fundações das estru- 
turas são actuadas por impulsos de aceleração 
a (t) em todas as direcções, podendo o conjunto 
de impulsos ser representado por um acelero- 
grama com várias componentes (fig. 1). 
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Fig. 1 — Acelerograma do sismo de Olympia, 
Washington, em Abril de 1965 


As condições de equilíbrio em cada elemento 
da estrutura são expressas por 


E+E+E=0O (1) 


di 4] 
| , | | | | A 
DON EM | Um Nm MAIA qua tade o proa É 


por ARTUR RAVARA * 


MADALENA QUIRINO ** 
MARIA TERESA PEREIRA DA SILVA *** 


SYNOPSIS 


In the studies on earthquake engineering that have been 
carried out at the Laboratório Nacional de Engenharia 
Civil for about 15 years, different techniques have been 
applied for the analysis of single and multidegree of jree- 
dom oscillators, linear and non-linear, 

The present paper describes the main characteristics of 
the systems thut have been considered, relating them to 
different types of structural behaviour, The formulation of 
the differential equations of motion as mell as the step-by- 
-step integration techniques are briefly reported. Examples 
of dynamic analysis of reinforced concrete structures acted 
by strong earthquakes are presented, 


sendo F,,F, e F, respectivamente as forças de 
inércia, de amortecimento e de ligação em cada 
elemento. 

Para oscilação segundo uma direcção xx e admi- 
tindo que o amortecimento é do tipo viscoso, caso 
usual em estruturas, as eg. (1) tomam a forma 


mx + Ce de ER) = 0 (2) 


se a massa da estrutura se puder considerar con- 
centrada num ponto (oscilador de um grau de 
liberdade), ou desdobram-se num sistema de 
equações 


 Xy) + 
Xx,) =0 (3) 


tt t 4 
mm, % CO Rgno e + 
e EM om » > 


jl=1,2,...,n 


se a massa da estrutura se considera concentrada 
em n pontos (oscilador de n graus de liberdade), 


(1) Comunicação apresentada ao Simpósio sobre as Teorias da Informação e dos Sistemas, Lisboa, Set. 1970. 
* Engenheiro Civil, Especialista do Serviço de Edifícios e Pontes, Chefe da Divisão de Dinâmica Aplicada. 
** Matemática, Especialista do Serviço de Edifícios e Pontes. 
“** Matemática, Estagiária para Especialista do Serviço de Edifícios e Pontes. 
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que no caso de edifícios correspondem geralmente 
aos níveis dos pisos. 

Durante um sismo o comportamento da estru- 
tura é determinado pelo seu movimento, relati- 
vamente ao movimento do solo. Desdobrando o 
deslocamento x em deslocamento relativo x, e 
deslocamento do solo x, 


x=x, —+x, (4) 


as eg. 2) e 3) tomam respectivamente a forma 


mx, + Cx, + F(x,)= — mx, (5) 
e 

m, Xi FC, (ksa o a RENA Xu) + 

+ E Quo Rara Covo 4 Rg) PR ER, X. (6) 


= 1,2 evse 
admitindo que o amortecimento depende apenas 
do deslocamento relativo. 


Designando a aceleração do terreno Xe por a (t) 
e suprimindo o índice r, dado que os desloca- 
mentos, velocidades e acelerações da estrutura 
que figuram nas eg. 5) e 6) são sempre relativos, 
obtém-se finalmente 


mx+ Cx+EG)=—m alt) (7) 


X) + 


sXn)=— m,a(t) (8) 


tr / t 
m,x,; + C, (XX ce. 
E. Mango » 


i=1,2,...,m 


A integração das equações 7) ou 8) permite 
obter os deslocamentos, velocidades e acelerações 
nos elementos da estrutura devidos à actuação 
de um sismo. Trata-se em todos os casos de equa- 
ções diferenciais ordinárias de 2.2 ordem, linea- 
res ou não conforme o tipo de funções F (x) con- 
siderado ; adopta-se geralmente para o termo 
de amortecimento expressões lineares da forma 


í t 
E DES numca 


No parágrafo seguinte referem-se os tipos de 
osciladores não lineares que têm sido estudados. 
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2 — OSCILADORES LINEARES E NÃO LINEA- 
RES 


Sejam dois osciladores de um e de n graus de 


liberdade, supostos com amortecimento viscoso, 
esquematizados na fig. 2. 


E 
| 


Fig. 2 — Osciladores de um e de n graus de liberdade 


+ n 


E 


E 


F, 


No caso das forças de ligação F (x) serem 
proporcionais aos deslocamentos relativos x, será, 
no caso de osciladores de n graus de liberdade 


n 


E => kz, (10) 


j=1 


equação que se reduz a F == Kx para n = 1. 
As eg. (7) e (8) tomam então respectivamente 
a forma 


mx+Cx+kKkx=—m a(t (11) 


m, x + 3 C, x, PA K, X; peça) 11 a (t) (12) 


j= j=1 


ou, no caso das eg. (12), adoptando a represen- 
tação matricial 


[m] (x) + [C] e) + [K] be] =— fm] a (t) 
(13) 
m, o| 


"a matriz de 
massa, 


0 m 


| 
sendo: [m] = 
L k 


[C] = [C,] a matriz de amortecimento e 
[K] = [K,|] a matriz rigidez, 
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designando por C, e K, respectivamente a força 
de amortecimento ou a força de ligação que se 
desenvolvem na massa i quando se aplica um 
deslocamento unitário na massa j e deslocamentos 
nulos nas restantes massas, ou seja, para x, = 
= 1 € X%ypy = O. 

Como é bem sabido, o comportamento dos 
materiais de que são feitas as estruturas de 
edifícios, nomeadamente aço e betão armado, só 
pode ser considerado aproximadamente linear 
até um certo nível de esforços ou de deformações 
aplicadas, como indica a fig. 3, em que F. e x, 
são respectivamente a força e o deslocamento 
elástico limite, acima dos quais o material entra 
em cedência. 


Xe Deslocamento, à 


se Deslocoment , 


a) Aço macio b) Betão simples ou betão armado 


Fig. 3 — Diagramas forças — deslocamentos 


No caso de sismos de pequena intensidade, a 
esquematização linear é válida e a análise dinâmica 
das estruturas corresponde à integração das eg. 
(11) e (12). Para sismos de forte intensidade 
porém, os numerosos estudos analíticos e experi- 
mentais realizados até à data [1] bem como a 
observação do comportamento dos edifícios em 
zonas onde ocorreram sismos intensos [2] indicam, 
de modo irrefutável, que os materiais que actual- 
mente constituem as estruturas são sujeitos a 
esforços e deformações que excedem largamente 
os limites elásticos, exibindo acentuada não 
linearidade acompanhada de efeitos de histerese 
mais ou menos pronunciados. 

No período mais intenso do sismo, de duração 
geralmente inferior a 30 segundos, os elementos 
estruturais são solicitados por ciclos de carga 
alternada de grande intensidade, cujos diagramas 
característicos, no caso do betão armado, têm o 
andamento indicado na fig. 4 [3,4]. 

O diagrama patenteia, além do decréscimo de 
rigidez quando aumentam as forças e os deslo- 
camentos, efeitos de histerese que aliás contri- 
buem largamente para a estabilidade da estru- 
tura. De facto, a área interior às curvas da fig. 4 
é proporcional à energia absorvida pela estru- 
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tura, a qual, adicionada à energia que corres- 
ponde às forças de amortecimento, deve equili- 
brar a energia transmitida pelo sismo à estru- 
tura. 


As: S8sqm 


Aç:158 sq m | 
E; :3,800ps: — 


t, 48,000 psi 
', :48,700 ps: 


Shearing force V.kips 


Fig. 4 — Diagrama característico de ciclos de carga 
alternada de grande intensidade em elementos 
de betão armado 


o) Diagrama elasto-plástico 
sem histerese 


b) Diagrams etasto-pidstico 


cj Disgrama be-linegr diDbicvgroma etasto-piastico 


de rigidez decrescente 


F 


| 
| T 


/ 


F| 
E 
K 
pdf 
x = x 


e) Disgromos bi-imecres com tronsições parabólicas 


Fig. 5 — Diagramas não lineares que têm sido adoptados 
em estudos realizados no LNEC. 
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Nos programas que têm vindo a ser desenvol- 
vidos no LNEC, a não linearidade tem sido 
expressa pelos diagramas do tipo indicado na 
fig. 5 a) a e). 

O efeito de histerese, que é considerado expli- 
citamente nos diagramas 5b) e 5e), pode ainda 
ser tomado em consideração quando se utilizam 
diagramas do tipo 5a), desde que se aumente 
o amortecimento de modo a englobar a energia 
dissipada por histerese [5]. 

É corrente a utilização de um diagrama cons- 
tituído por curvas tipo Ramberg-Osgood [6] indi- 
cado na fig. 6. Berg [7] e Jennings [8], entre outros 
autores, têm-no utilizado. Na prática difere pouco 
dos diagramas 5e), os quais têm sido preferidos 
nos estudos realizados no LNEC devido a conte- 
rem como caso particular o comportamento 
linear, que se verifica em grande parte das estru- 
turas, mesmo durante um sismo intenso. Com 
efeito a plastificação restringe-se a algumas zonas 
e como a análise se complica consideravelmente 
devido à não linearidade, a utilização dos diagra- 
mas 5e) conduz a uma economia global de cálculos 
importante, relativamente ao diagrama Ramberg- 
Osgood. 


x+xo F+ F+4 Fa” 
ido À à EV => - 
2x, 2F, 28 4 


« f A=A, TO (EE 


2x, 2Fs 2F 


Fig. 6 — Diagrama tipo Ramberg-Osgood 


Um outro tipo de não-linearidade que pode 
ser necessário considerar é o devido a efeitos 
gravíticos, quando os osciladores sofrem grandes 
deslocamentos [9.] 

No caso do sistema considerado na fig. 7, 
o efeito gravítico pode traduzir-se aproximada- 
mente por uma redução no termo de restituição 
da eg. 7), vindo então 


nen de EA ho FO) 1 |=-mato (14) 


sendo P o peso da massa m e 1 a respectiva 
cota. 
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o (t) 


Fig. 7 — Efeito gravítico num oscilador de um grau 
de liberdade 


3 — MÉTODOS DE CÁLCULO 


Na integração das equações que descrevem 
o movimento dos osciladores utilizaram-se dois 
métodos de Runge-Kutta de 4.2 ordem [10, 11), 
aplicáveis a sistemas de equações diferenciais 
ordinárias de 1.2 e de 2.2 ordem. 

A integração é feita por um processo «passo 
a passo», sendo a aceleração a (t) dada para um 
conjunto discreto de valores da variável t. A par- 
tir das condições iniciais do movimento e dos 
valores da aceleração, são sucessivamente calcula- 
dos os acréscimos das grandezas que intervêm 
nas equações do movimento, para intervalos de 
tempo Át. 

No quadro 1 apresentam-se as principais carac- 
terísticas dos métodos utilizados. 

Qualquer dos métodos se revelou satisfatório, 
no que diz respeito a precisão e rapidez de con- 
vergência. 

Em sistemas elasto-plásticos e bilineares a inte- 
gração é desdobrada nos domínios corresponden- 
tes aos ramos das rectas dos diagramas 5a) a d). 
As condições de transição são determinadas pelo 
valor do deslocamento e pelo sinal que toma 
a velocidade. 

Nos sistemas bilineares com transições para- 
bólicas sucede o mesmo, com a diferença de que 
a rigidez, que naqueles sistemas é constante den- 
tro de cada domínio, varia nestes de ponto para 
ponto ao longo dos ramos parabólicos. 

O acelerograma é em geral definido para 
valores do tempo espaçados de 0,05 seg; no caso 
de serem necessários mais valores obtêm-se por 
interpolação linear daqueles; o passo de integra- 
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QUADRO 1 


Características dos métodos de Runge-Kutta utilizados 


Desi- 
gnação Runge-Kutta Runge-Kutta-Merson 
Ordem 4.a 4a 
Campo | sistemas de equações diferençiais ordiná-| sistemas de equações diferenciais ordinárias 
de rias de 2.º ordem de 1.2 ordem 
aplicação 
sistema na forma matricial: Passa-se do sistema: 
X,=5S,—CX,—K,DX q , : 
X,=F(X,X,) x,;=g(x, dd n 
Yaoi 
Formu- EEN ES SUNT E ai 
lação a por: t=X;X,= Ya; p1/X; dx Yoi+2 
para o sistema: 
“ ! dy i 
DXa)= (x); bx j=1x) ea Mme Faia (Ya a) 
j ni: 
[Dx|=[4x] dYscsã 


dx Foto (Yoiro Y2iy10Y9i-10Y2i-s) 
= 1,2 «sa; À 


DX=At. 6, Yi =Yia Hky,4k, + ko);/2 
Cálculo | X,=X,, +4Át.6, ne 1,2. cos À 
sendo X,, o valor de X anterior a X, sendo y;, O valor de y anterior a y, 
At h 
é ui Go ud ku H04Yai) 
do,=(k,+2k,;4+2k;+k,... 4k,) 6 h 


h 
mato aco | 


em que: 
k,=F(t,X,X,) h k 
o 2H) El o) a) | 
At At At >i Ep = pd E da 2 
Algo- ko=E(t+ Pica E E de 7 a 
ritmo de . k = Deja do(yas 8) | 
cálculo ks=F(t+ At mpLy 4 (A t)" k,, ni ai 2 Ya 8 
2 2 4 
BE. E el 3 
+ 2 a) a ie (ot? E tok), | 
2 
k=F(* FALX+HAEX + (At) ho; em que h==At (passo de integração) 
4 
X+At. ks) 
| Não calcula o erro. O passo de integração | O erro E, é dado por 
Precisão i 
E deve ser suficientemente pequeno para 9k, 5 
esa garantir boa aproximação. =|| k— 2. +4k,— 
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ção depende das características dinâmicas dos 
osciladores e a sua ordem de grandeza é a do 
período próprio mais pequeno. É assim corrente, 
em sistemas de 10 a 15 graus de liberdade, cor- 
respondentes a edifícios com este número de 
pisos, um passo de integração de 0,005 seg. 
O passo de saída dos resultados pode ser qual- 
quer múltiplo do passo de integração. 

A descrição pormenorizada das técnicas apli- 
cadas consta de diversos relatórios [12 a 15]. 


4. APLICAÇÕES AO ESTUDO SÍSMICO DE 
EDIFÍCIOS 


4.1 — Sistemas de um grau de liberdade 


São aproximadamente redutíveis a osciladores 
de um grau de liberdade estruturas de edifícios 
dos tipos indicados na fig, 8 a) e b) em que a 
deformabilidade do piso ou dos pisos inferiores 
é muito superior à dos restantes pisos. 


Pisos que 
permanecem 


intoctos é 


Pisos em que 
hó cedência 


a) b] c) 


Fig. 8 — Estruturas redutíveis a osciladores de um 
grau de liberdade 


Transformam-se ainda em sistemas deste tipo 
edifícios de vários pisos com estrutura aporticada, 
preenchida ou não por alvenaria, ao fim de certo 


cm 


— a) 
4 -— cs 


ns, 


t 
f 
| 
Teperiodo | 
n-omortecimento relativo | 


io «u aa é id A 


d 
“f 
u 


Deslocamen: 


tempo de actuação de um sismo intenso. Com 
efeito, a cedência restringe-se em geral aos pisos 
inferiores, que adquirem assim uma deformabi- 
lidade muito superior à dos pisos restantes, 
fig. 8 c). 

Por outro lado uma das técnicas mais potentes 
para a análise dinâmica de sistemas de vários graus 
de liberdade em regime linear é a análise modal 
[16], na qual se determina a resposta do sistema 
por sobreposição da resposta de tantos osciladores 
de um grau de liberdade quantos os modos de 
vibração que contribuem para a resposta do 
sistema global (!). 

Na análise de osciladores de um grau de 
liberdade têm sido utilizados programas de cálculo 
dos seguintes tipos (!!), considerando-se em todos 
os casos amortecimento viscoso : 

— Oscilador linear: o estudo de sistemas deste 
tipo reveste-se de grande interesse, devido a ser 
usual caracterizar um sistema por meio dos 
valores máximos do deslocamento, velocidade e 
aceleração que se desenvolvem em osciladores 
lineares com amortecimento viscoso. Designam-se 
por espectros as relações entre aqueles valores 
máximos e as características dinâmicas dos osci- 
ladores. 

Na figura 9 apresentam-se os diagramas de 
resposta de osciladores do tipo referido, actuados 
pela componente N-S do acelerograma do sismo 
de El Centro, California, de 18/5/1940. 


(1) É corrente considerar os modos correspondentes a 
frequências próprias inferiores a 5 Hz. 

(11) O Centro de Cálculo do LNEC está presentemente 
equipado com um computador NCR Elliott 4130. 


| 
| 
mo ; 
(é E 10 
Tempo, S£9 


Fig. 9 — Deslocamentos de osciladores lineares de um grau dé liberdade actuados pela componente 
N-S do sismo de El Centro de 1940 
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Velocidades (mms) 


Na fig. 10 apresentam-se os espectros de des- 
locamento, velocidade e aceleração para a com- 
ponente E — O do sismo de Lisboa de 28/2/1969. 


200 


BE SRS SD A É 


x o 
“ah ExENTX 
SL PR NV 


X| 7 “e 


Rest PR DES 
AS, 68 


frequência, Hz 


Fig. 10 — Espectros de deslocamento, velocidade e 
aceleração da componente EO do sismo de Lisboa 
de 28/2/69 


— Oscilador bilinear; prepararam-se três pro- 
gramas, destinados respectivamente à análise 
de osciladores bilineares (contendo como caso 
particular o oscilador elasto-plástico), de oscila- 
dores bilineares com rigidez decrescente (fig. 5d) 
e de osciladores bilineares tendo em conta o 
efeito gravítico, com o comportamento descrito 
pela eg. 14). Os programas têm sido largamente 


aplicados no estudo de diversos tipos de estru- 


=] 


Deslocamento relativo x/x, 
Q 


E 
“3 
-4 
jr 
ua 
O 1 2 3 4 5 


— Oscilador brlinear. w,= 2H72;17:0,05,x,5:2,5cm;r:0,12 
4[--— Oscilador elasto-plástico:w = 2H2,1:0,05,x9=:2,5cm;r=O0 


turas. As características dos diagramas estão 
directamente relacionadas com as propriedades 
dos materiais e com a própria constituição das 
estruturas. Assim, por exemplo, o diagrama elas- 
to-plástico é adequado para esquematizar o 
comportamento de estruturas metálicas e de 
estruturas de betão armado com pequena per- 
centagem de armadura (4). Ao aumento de per- 
centagem de armadura corresponde o aumento 
de inclinação do ramo II do diagrama bilinear. 

O efeito do esforço axial sobre os pilares das 
estruturas também se traduz, dentro de certos 
limites, por uma variação de inclinação do ramo 
Il do diagrama (4). 

Na figura 11 compara-se a resposta de um 
oscilador bilinear com a de um oscilador elasto- 
plástico, actuados por um sismo intenso. 

— Oscilador bilinear com transições parabólicas : 
dispõe-se de programas de cálculo para os dia- 
gramas indicados na fig. 5 e) e f). Na maioria 
das aplicações os dois diagramas são prática- 


“mente equivalentes, correspondendo, no caso de 


estruturas de betão armado, a um aperfeiçoamento 
considerável relativamente à esquematização bi- 
linear. Os programas têm sido aplicados no estudo 
do comportamento de um grande número de 
osciladores, com o principal objectivo de rela- 
cionar as principais propriedades dos materiais e 
das estruturas com as características da res- 
posta (4). A fig. 12 exemplifica o tipo de infor- 
mação colhida. 

Todos os programas referidos neste parágrafo 
fornecem os valores dos deslocamentos, veloci- 
dades e forças de restituição correspondentes 
a valores do tempo dados. No caso de osci- 
ladores não lineares, são assinaladas todas as 
transições dos diagramas de restituição. 


CARACTERISTICA BI-LINEAR 


we = Frequência própria 


n = Amortecimento relativo 


. t9B a /A9a 
é “ga e m 


Tempo ,Sseg. 


Fig. 11 — Deslocamentos de um oscilador bilinear e de um oscilador elasto-plástico actuados pela componente 
N-S do sismo de El Centro, de 1940 
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Deslocamentos máximos cm 


—— — uu 


Pilares: 
B 300, A4O 

G = 60 Kgt em? 
b=40cm,h,:60cm” 


— ) mãx atingido 


— — D max, admissivel 
(Ep: 3,5%he VE; =10%00) 

e - Despreza-se o e axial 

O - Considera-se a influência | 
do e axial | 


10 1,5 2,0 
Percentagem de armadura TÚ, 


Fig. 12 — Resposta de estruturas redutíveis a osciladores bilineares com transições parabólicas. 


4.2 — Sistemas com n graus de liberdade 


Prepararam-se dois programas de integração de 
sistemas de equações relativos a estruturas de 
dois tipos muito diferentes: 


Estruturas aporticadas de pisos rígidos 


Situam-se nesta categoria as estruturas em que 
a rigidez dos pisos é muito superior à dos pila- 
res, que se podem considerar, quando a estru- 
tura é actuada por forças horizontais, aproxima- 
damente encastrados nos pisos. 

A força de ligação em cada massa m, depende 
então apenas do deslocamento das massas con- 
tíguas m, , em,,, pelo que o sistema (8) toma 
a forma 


mx +2nw mx + EFi(x—x |) + 


A. Fº ,—x,.))=—m,a(t) (15) 


e TZ orecy NR 


sendo F, (x,—x;,. )) e E; (x,—x;.1) as forças 
de ligação respectivamente à massa inferior e 
superior a m, e adoptando um tipo muito sim- 
ples de amortecimento viscoso, com um valor 
relativo constante 7 e proporcional à frequência 


própria fundamental w.. 
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O programa de cálculo para a integração deste 
tipo de sistemas (14) considera relações forças- 
-deslocamentos de tipo misto, bilinear com tran- 
sições parabólicas, indicado na fig. 5 e). O pro- 
grama opera a partir das características da estru- 
tura e das acelerações sísmicas. Dos resultados 
constam, com a disposição indicada, os seguintes 
valores : 


tempo, t, 
E neRe Velocidade Forças de lidação 
x x F: (x,) pe (x ) 
| 1 pj 11X, —Xga 
x x E ) E ) 
9 9 2 Mo —Xj 9 Mg — Xa 
1 
Xn Xn Fa — Hg =2) 0 


Os valores de t, são, além dos corresponden- 
tes a igual espaçamento da resposta, todos aque- 
les para os quais se verifiquem mudanças de 
ramo na característica forças-deslocamentos em 
qualquer dos pisos. 

Apresenta-se na fig. 13 um extracto da resposta 
de um oscilador de 8 graus de liberdade (estru- 
tura de edifício com 8 pisos), com os resultados 
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3.200 


0.37 
0.68 
0.93 
1.16 
1.37 
1.55 
1.89 
2.09 


5.97 
11.06 
16.37 
21.85 
26.29 
29.25 
32.76 
35.83 


12.49 
2.73 
7.36 
5.56 
4.48 
3.75 
2.88 
1.55 


expressos em segundos, centímetros e toneladas 
força. 


a É, 
— 7.36 
= 3.96 
— 4,48 
=3,75 
— 2.88 
=1.55 

0.00 


Fig. 13 — Extracto da resposta de um oscilador de 8 
graus de liberdade 


O programa permite a análise de osciladores 
com um número máximo de graus de liberdade 
igual a 100. Tem sido utilizado em estudos de 
indole semelhante aos referidos no parágrafo 
anterior, apresentando-se a título de exemplo, 
na fig. 14, resultados do estudo sísmico de um 
edifício de 8 pisos. 


Estruturas de tipo consola ou misto 


A evolução da tecnologia da construção de edi- 
fícios altos de betão armado está conduzindo 
a soluções estruturais cujo comportamento se 
afasta do das estruturas de pisos rígidos referi- 
dos no parágrafo anterior [2]. Com efeito a ten- 
dência geral é reduzir a rigidez dos pisos relati- 
vamente à dos elementos verticais, pilares, pare- 
des e núcleos resistentes. Daqui resulta que 
a deformabilidade sob a acção de forças hori- 
zontais se pode decompor em flexão, corte e rota- 
ção da base, como indica a fig. 15 [17]. É então 
sempre possível, dentro de certas hipóteses simpli- 
ficativas correntes em engenharia sísmica [18], 
definir uma consola com características geomé- 
tricas e mecânicas que conduzem a uma defor- 
mabilidade igual à da estrutura inicial. As equa- 
ções (13) são adequadas para descrever o com- 
portamento desta consola equivalente, com mas- 
sas concentradas em n pontos, actuada por uma 
aceleração na base, sendo a matriz rigidez K dada 
por 
E ed bes bd As be=b” Eb (16) 


8 To CARACTERÍSTICAS DA ESTRUTURA 

a | — Tensão normal permanente na base des pilares 1 — 2: 56 kgf cm dida 
| — Pilares: 

Ss 1—2:40>40 cm;.=0,75", 5—6:30>X40cm;n,=0,5 “q 
+ nje 5 2— 3:40>(40 cm;w =9,75")p 6—7:30x40cm;o.=0,5 "4 
É A 3—4:40><40 m;w =0,75“h 7—8:30><30cm;w.=0,5 º/, 
c 4— 5:30>40 cm;u.=0,5 “py 8B—9:30> 30 cm;u. = 0,25 "4 
81 ” — Frequência própria : 1,40 HZ 

e 2 -— Amortecimento viscoso: 5 “/, 

nO + 


FORÇAS MAX ,tf 


2,27 
—— 
2as 
———s 
1,38 
fcáno 7 
139 COEF SISMICO=022 
3.60 
3 
1,68 
4 
1,96 
——e 


DESLOCAMENTOS MAX., em 


—— MA rimo atingido 


---- Múximo admissivel, correspondente o 
extensões no betão e no aço de 35 ue 
10 es 


—— Elostico limite 


Es 19 12em 


Fig. 14 — Forças e deslocamentos máximos numa estrutura de pisos rígidos 
sob a acção de um sismo de forte intensidade 
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-— 


DESLOCAMENTO FLEXÃO CORTE pOTAÇA 
Da sat 
TOTAL SAS 


Fig. 13 — Deformabilidade de uma estrutura mista 


sendo 4; a matriz de deformabilidade, correspon- 
dente à consola encastrada, b a matriz que rela- 
ciona forças horizontais com momentos flectores 
nos pisos, f a matriz de deformabilidade dos 
tramos que “constituem a consola equivalente, 
relacionando momentos flectores e rotações nas 
extremidades e As a matriz deformabilidade cor- 
respondente à rotação da base (18). 


A matriz f por sua vez é da forma 


-t 
» | 
f= , | (17) 


= 


em que 


[£,] ia [£]m E [£,] (18) 


referindo-se os índices M e T respectivamente a 
deformabilidade por flexão e por corte. As rela- 
ções entre f, e as características geométricas e 
mecânicas de cada tramo são bem conhecidas [19]. 

Preparou-se um programa designado por ON- 
LINGRAL [15] para a análise dinâmica das con- 
solas equivalentes, considerando o comporta- 
mento não linear correspondente ao diagrama 5f) 
para a deformabilidade por corte e por flexão, 
independentemente uma da outra. 

O programa efectua o cálculo dos desloca- 
mentos, velocidades, forças de ligação, momen- 
tos flectores e esforços transversos ao nível a 
que estão concentradas as massas, apresentando 
os resultados sob uma forma semelhante à indi- 
cada na fig. 13. Está em curso a análise do compor- 
tamento de uma gama de edifícios de 8 a 20 pisos, 
para a determinação da influência das caracterís- 
ticas estruturais na resposta. As características que 
estão a ser consideradas são, além da geometria 
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global do edifício, a frequência própria, o amor- 
tecimento, a configuração dos modos de vibração 
e muito especialmente os tipos de não linearidade 
devidos aos elementos estruturais e aos elemen- 
tos de enchimento de alvenaria. 

A fig. 16 apresenta resultados obtidos na aná- 
lise de um sistema de 8 graus de liberdade, cor- 
respondente a uma estrutura aporticada de betão 
armado preenchida por elementos de alvenaria. 


Iser “sec S sec E 
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| mm DENDMO 
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mto eee em e em 


tm 


| 4º" STORY 
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Fig. 16 — Resposta de um oscilador de 8 graus de 
liberdade [17] 


Qualquer dos programas descritos em 4.2 con- 
tém como caso particular os sistemas lineares. 
O último dos programas referidos é o mais geral 
permitindo a análise de quase todos os tipos 
de sistemas considerados ao longo do trabalho, 
para diagramas de restituição de tipo misto. Aliás 
a alteração do tipo de diagrama característico 
implica apenas pequenos ajustamentos no pro- 
grama. 


5 — CONCLUSÕES 


A Engenharia Sísmica constitui um vasto campo 
de aplicação de técnicas numéricas para a inte- 
gração de equações e de sistemas de equações 
diferenciais ordinárias de 1.º e 2.2 ordem, linea- 
res e não lineares, correspondentes aos modelos 
matemáticos de diversos tipos de estruturas. 
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Ao longo de cerca de 15 anos de actividade 
naquele domínio foi elaborado no LNEC um con- 
junto de programas de cálculo destinado a cobrir 
os principais tipos de estruturas de edifícios. 
O planeamento e a realização desta actividade 
constitui um trabalho de equipa e exige um 
suporte experimental importante, que lhe é pro- 
porcionado pela realização de estudos sobre 
modelos e pela observação do comportamento 
das construções sujeitas à acção de sismos de 
grande intensidade. 
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Indústria Automóvel. 
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OBJECTIVO MAIOR. 
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APLICAÇÕES DA GEOESTATÍSTICA (*) 


RESUMO 


Em três problemas práticos, pretende-se provar que os 
metodos da qgeoestatística, na medida em que conferem ao 
depósito mineral wna certa «estrutura» que lhe é dada pelo 
mudo particular como ocorreu a mineralização, conduzem a 
menores erros do que um tratamento pela estatistica clássica, 

Esses problemas são: 

1. Definição da superficie e volume mineralizado, Erro 
deste cálculo, 

2. Variância de extensão de uma amostra à sua zona 
de influência, 

3. A Krigeagem — ou seja a procura do melhor estima- 
dor para um certo parâmetro (aquele que minimize a 
variância definida em 2,), 


1. DEFINIÇÃO DA SUPERFÍCIE DE UM 
JAZIGO 


No reconhecimento de um depósito mineral e 
no cômputo das suas reservas, pretende-se em 
primeira análise defini-lo em planta — estimar a 
sua superfície e calcular o erro dessa estimativa. 
Temos uma rede regular de sondagens que nos 
dão uma informação —-a sondagem encontra a 
formação ou não. Define-se uma função k (x) — 
chamada indicatriz — tal que toma o valor 1 se a 
sondagem encontra a formação e o valor O se 
não encontra. 

Formemos agora o produto k (x). k(x + h), 
capaz de nos dar já uma ideia da estrutura do 
jazigo (sendo x uma coordenada de espaço, esse 
produto é o produto dos valores da indicatriz em 
dois pontos distantes do vector h). Então k (x). 
k(x + h) vale 1 se os dois pontos pertencem 
ambos à formação e vale o no caso contrário. 

Integrando em x (*), obtemos o chamado cova- 


(*) — Em tudo o que se segue, x representa um vector 
a n coordenadas. Pode ser uma coordenada linear, de 
superfície ou de volume e o integral correspondente pode 
ser duplo ou triplo. 


por HENRIQUE GARCIA PEREIRA 


Engenheiro Quimico-Industrial 


SYNOPSIS 


In three practical problems, it is intented to prove 
how geostatistic methods, as they take in accout the 
«structure» given to the ore deposit by its particular mine- 
ralization, lead to better sigmficance levels (to amaller 
errors) than classical statistics. Those problems are: 

1) Definition of the deposit surface and volume. Error 
of this calculation, 

2) Extension variance of a sample to its influence 
zone, 

3) The Nrigeage — the search of the best estimatron of 
a certain paramether minimizing the variance treated 
in 2.). 


riograma geométrico associado à estimação, que é 
uma função característica de um conjunto de 
métodos da gecestatística denominados métodos 
transitivos que são adequados à resolução de 
problemas de transição, de fronteira: 


-+ OS 


fi) es | Ok te + h) dx 


« 
— OO 


Se dx representar um elemento de superfície, 
k (h) é a Medida da superfície intersepção de S 
com 5 transladada de -h (a integranda só existe 
dentro desse domínio e aí vale 1): 


k (h)=Med SNS, 


Vejamos algumas propriedades do covario- 
grama: 


a) K(h) =K(—h) 
b)K(o) =sS 


o) K(h) =K(o) 


(*) Texto base de um colóquio realizado no Laboratório de Preparação de Minérios do 1.5.T. e integrado no Semi- 
nário de Geoestatística e Automação do Núcleo de Tratamento de Minérios do C.E.E.N. 
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Estas propriedades estão patentes na interpre- 
tação gráfica da Fig. 1: 

A função é simétrica, parte da origem com o 
valor $ e anula-se para um valor a, dito amplitude 
(«portée») que coincide com o diâmetro aparente 
da superfície na direcção « do vector h. Vemos 
portanto que existe um covariograma (no caso 
geral da superfície ser irregular) para cada 
direcção do vector h, veículo da translaçção. 


po 


S 


Fig. 1 


Outra propriedade importante do covariograma 
é o comportamento na vizinhança da origem — se 
h é suficientemente pequeno, o covariograma é 
linear : 


K (h)=K(o0)—D,h 


onde D, é o diâmetro aparente da superfície 
perpendicularmente a «x, e portantoD,h é a 
superfície infinitesimal que é necessário retirar a 
S (que é K (0)) para obter SS". 


K (h) — K 
Então, Dia h da 


ou seja, o diâmetro D, é a derivada do covario- 
grama na origem e este pode obter-se de D, 
por integração. 
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Utilizando o covariograma geométrico como 
ferramenta, formulemos o problema que nos 
propusémos : 

Consideremos uma formação mineralizada 5 
reconhecida através de sondagens dispostas 
segundo uma malha regular rectangular de lados 
a, e a,. Um estimador dessa superfície será (*) : 


x E 
PRP Rita Ys Pag) 
| Pp 


onde K(x, y) (com x=x,+ajie y=y,+ 
“+-a,p), é a indicatriz, agora explicitada a duas 
dimensões, que vale 1 se a sondagem centrada 
em (x, y) encontra a formação e o no caso con- 
trário. 


Portanto, o duplo somatório dá o número n de 
sondagens positivas usando esta rede com a ori- 
gem aleatóriamente implantada em (x, , y,). Assim, 
S será uma variável aleatória (por força da impla- 
tação aleatória da origem — ao iniciarmos uma 
campanha de sondagens não temos nenhuma 
ideia de onde centrar a origem da rede e o valor 
obtido para o estimador S depende do ponto 
escolhido para origem). O estimador S é «bom» 
se a esperança matemática E(S)=S e o erro 
cometido é S — S. Vamos tratar o problema a uma 
dimensão. Para isso suporemos que o erro come- 
tido na estimação de uma linha por uma série 
de sondagens descontínuas é independente do 


erro que se comete estendendo o valor obtido 
anteriormente à sua zona de influência. 


(*) X representa sempre o valor estimado da variável 
cujo verdadeiro valor (desconhecido) é X. 
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Então um estimador do comprimento da linha i é 
L= a, >, K (x, +Ha i, yo + aopP) 
p 
S pode escrever-se : Ss à > L 
e S — S =5-—a, bz L, 


somando e subtraindo a, b> L, , obtemos 


S-S=(S-a, 2 L) + a E (ut) (1) 


O primeiro termo, S — a, o L,, é o erro de 


fatias (erro que se comete admitindo que já se 
conhece o verdadeiro valor dos L, e que traba- 
lhamos a uma dimensão — ao longo dos p — como 
se cada comprimento L, estivesse acumulado num 


ponto) eo L, — L, é o termo de linhas (erro que 
se comete por termos amostrado o comprimento 
L;, por sondagens descontinuas). 

Assim, o erro em S será a soma do erro de 
fatias com o produto de a, pelo erro de linhas. 
Vamos calcular cada um destes erros independen- 
temente. 

O erro de linhas pode obter-se pelo seguinte 
raciocínio : 

Suponhamos um jazigo linear de compri- 
mento b, a estimar, reconhecido por uma malha 
regular de sondagens de passo a. A implantação 
aleatória da origem leva a que o jazigo possa 
admitir n sondagens positivas e então b = na 
ou n-t 1 sondagens positivas — e então b = 
(n + 1) a. Se for Ea o comprimento que o jazigo 
tem a mais do que n« (com : << 1), o verdadeiro 
valor b (desconhecido) é b == (na + E ajea pro- 
babilidade de implantação da origem entre O e 
(1—:) a é (1 —:). A probabilidade de implanta- 
ção entre (1I—:)a ea és. Á primeira possibili- 
dade corresponde um valor na e à segunda, 
(n + 1)a. 

Calculemos a esperança e a variância do esti- 
mador : 


E (b=(n+1)a.ti na(l-s)=na+ta=b 

22 0) es Eb be E tp = E (|? 

52 (b) =(n+1) atietnal. (1-4) — 
. — (na + ta) 

cb) =(:— ea 
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Xg => bzna 
| Xxg => b=(nel)a 


E a 


(nella 
0 a 2a na 


/L// LLC LL LL LL LL LL LL 
(1-Eja 
GD aa 


Fig. 4 


O valor médio desta expressão, admitindo uma 
distribuição de probabilidade uniforme, é 
| 
sb) ea a“ [ (+) di == E a 
1 j (o) 
o 

Este valor que obtivemos por um raciocínio 
puramente probabilístico foi obtido por Mathe- 
ron” por tratamento sobre o covariograma. Con- 
cluiu que os termos em — r do covariograma dão 
para a variância termos em 1/6 a” (o que nos 
leva a concluir que o covariograma da estimativa 
de um segmento por uma rede de sondagens 
regular é linear, como é fácil de verificar), os 
termos r” log r dão termos em 0,0609 a” e deum 
modo geral, os termos em r* dão para a variância 
termos em A; r”'! (com À impar). A estas rela- 
ções que permitem passar directamente do cova- 
riograma para a variância, chamou Matheron 
princípios de correspondência. 

Quanto ao erro de fatias, esse não pode 
ser calculado sobre o covariograma geométrico 
associado à estimação (que é linear, da forma 
k (h) =k(o) — Dh ), mas sim sobre este depois 
de ter sofrido uma «montée» de primeira ordem 
de modo a ter em conta a acumulação das diver- 


sas linhas nos pontos distanciados de a,. Sobre 


R 2 : di 
a teoria da montée, ver!” e !”, Para as aplicações, 


o problema põe-se assim: Dado um teor em 
volume t (x, y, Z) — variável regionalizada dotada 
de um covariograma a três dimensões — se inte- 
grarmos os teores ao longo dos ZZ, obtemos um 
teor em superfície ou acumulação no ponto (x,y): 


z 


bo) 
- 


t (x, y) = | t (x,y, z) dz 
I 
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t (x,y) será uma nova variável regionalizada 
dotada de um novo covariograma que pode ser 
obtido do anterior por «montée». 

Matheron demonstra que a montée pode ser 
efectuada directamente sobre o covariograma 
(sem ser necessário integrar sobre a própria 
variável) e dá para esta operação uma série de 
princípios de correspondência (paralelos aos já 
escritos para as variâncias) : 

p=» qa 1 
logr — rr 
—r —r log r 


Então para o nosso problema, se o covario- 
grama é linear, o termo de linhas dá 1/6 D as 
e o termo de fatias obtém-se do covariograma 
monte — que é D hº log h — pelos princípios de 
correspondência escritos para a variância e dá 
Dx 0,0609 a; . 


Então o? (5)=D (0.0609 a, E ç a, 3) como 


+ 


se vê de [1] 


Fazendo D = N, a, =N, a, e dividindo por S* = 
= (n a, O obtém-se 


' 2 
ER gas N, + 0.0609 2) [2] 
G? nº 1 6 Ni 


fórmula importante para as aplicações, onde n é 
o número total de sondagens positivas, N, e N, 
são o número de sondagens máximo nas direc- 
ções p e i respectivamente (Fig. 3) com N,-N,. 

Sobre a aplicação concreta da fórmula [2] 
convém fazer as seguintes observações : 


1 — Se as direcções 2, e “2 de maior e menor 
dimensão do jazigo não coincidem com as direc- 
ções de propagação da malha, a fórmula [2] não 
é directamente aplicável e teremos de corrigir 
(segundo SERRA (3)) a anisotropia da malha em 
face do jazigo através de uma afinidade de eixo 
«, e módulo D,/D, . Então a malha (por hipótese 
rectangular) torna-se num paralelogramo de lados 
b e E e ângulo agudo 6 (Fig. 5) e a variância 
na estimativa da superfície é 


) 3924 |/ 
2=K, A“WVD,D,, 


aos 


onde ki é uma função de 9, & e E” dada em 
ábacos em (3) e A=-bb' é a área da malha 
antes de sofrer a transformação. 

Ee D; nação 


AÇ:—. . AB 
Dz 


Fig. 5 


No caso de malhas regulares, K, toma a forma 


simples 
, 1/92 fr 3 
K=1(, ) + 0.060 ( —) 
6 No senb 6 sen 


onde se reconhece a estrutura da fórmula [2] 
apresentada anteriormente, mas corrigida para 
anular o efeito da anisotropia. 


2 — No caso do jazigo conter pontos interiores 
de sondagens nulas, n continua a ser o número 
total de sondagens positivas, mas a contagem 
dos N, e N, faz-se como se segue: 

Considera-se 2D, como todo o perímetro mine- 
ralizado (que é um múltiplo da malha) tanto no 
interior como no exterior do jazigo paralelamente 
a uma das direcções de alongamento e 2D, para- 


lelo à outra. Então N, = a eN, = D,/a (se 
a 


a malha for quadrada e a direcção de propa- 
gação da malha coincidir com a dimensão máxima 
do jazigo). 

Apliquemos a noções apresentadas ao cálculo 
do volume mineralizado de um jazigo sedimentar 
de ferro, representado esquematicamente em 
planta na Fig. 6: 


Os números que afectam cada sondagem repre- 
sentam a altura mineralizada em metros e a 
malha é rectangular de 100 200 m. 

Calculemos primeiro um estimador da super- 
fície e o respectivo erro: 
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E; 
— a ———e 
| 200m , 
Fe 
+0.9 o 
Oo 
+ 0.4 +0.8 
Lá 
+1,2 +2.8 +2.0 
+0.5 EA +3,7 + 2,0 
é 
+ 1,4 + 1,2 + 0,5 +2.0 
+ 1.6 +0,7 +2.2 +0.5 +2,7 E, 
+ 1.7 
+ 0,2 +2.8 + 19 
+0.3 
Fig. 6 


Escolhidas as direcções «, e x, de maior e menor 
desenvolvimento do jazigo, teremos 


n==25 
a == 100m 2 D,=30 a 
x 
N,= 20X100. 15 
2 x 100 
b=200m  2D—-24b 
BE cam ER 
2 >< 200 
; Ea 1 [12 15º 
a fórmula (2) dá = as | e + 00605 =| 
” 2 LO 12 
q o: = 
— = 0.005 == 009 DO aa A 
9” = 5 


O estimador da superfície é 


* 
S = 25 (200 x 100) = 500 000 mº 


7 + . . A * 
— com um nível de significância de 
2 95 0 


dr 
4=5 2 
S = 500000 + 14% ou S= (500 +70) 10º m* 


Já conhecemos a área do jazigo e calculámos o 
erro dessa estimação através da fórmula [2], 
desconhecida da estatística clássica, que faz inter- 
vir as noções de covariograma e de «montée» 
e confere assim à variável em estudo um certo 
aspecto estrutural (não é uma variável puramente 
aleatória, é uma variável regionalizada). De resto 
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notemos a dificuldade que existe em dar um tra- 
tamento probabilístico a uma variável do tipo 
tudo ou nada — a sondagem encontra ou não a 
formação — e este acontecimento não pode vir 
afectado por nenhuma probabilidade. 

Para calcular o volume do jazigo vamos dar 
à variável «altura mineralizada» um tratamento 
diferente. Se for h, a altura de minério obtida 
em cada sondagem, uma estimativa do volume 
mineralizado será 

Vs 2 h, 


onde s é a área de um pequeno rectângulo de 
área ab, centrado no ponto x. 


onde n é o número de sondagens. 
O volume é pois igual ao produto da área 
total pela altura média. 


Supondo S e h duas variáveis independentes, 
2 ad 
v s h 
= ta 


Já calculámos o pelos métodos da geoestatis- 
tica. Vamos calcular c pelos métodos da estatis- 


tica clássica. Vamos considerar assim os valores 
das espessuras mineralizadas como a realização 
de uma certa variável aleatória, com uma certa 
lei de distribuição de probabilidade constante 
(mas desconhecida) : 


=-— =1.52 e Ehi=g7405 
n Í 
ZW o 
62, = | o hó= 0.64 
n 
[r) 0,64 
9 hi 
q. = O Ca SE Um 
d E 0,0256 
o. 0.0256 
— == (0,011 
hº 2.32 
O 2 É 
—— == 0,105 e —— = 0.21 
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ai DAM + pa 
h=h+ 
h 
ou h = (1.52 + 0.32) m 


== 1,52 + 21h 


Como vemos, o erro relativo na altura (com o 
mesmo nível de significância de 95"/9) é de 21º/o, 
enquanto que o erro na superfície é apenas 
14,0º/. Esse erro é inferior na estimativa da 
superfície S (apesar das irregularidades do jazigo) 
porque conferimos a essa variável uma certa 
estrutura (ao aplicarmos os métodos da geoestatis- 
tica). E essa estrutura é real, dado que o jazigo 
foi formado de uma certa maneira (explicada pela 
geologia) que não pode ser considerada aleatória. 
O cálculo da variância pela estatística clássica 
(erro do valor médio das espessuras) conduziu a 
um estimador «pior ». 

Vejamos o peso dos dois erros no erro volume: 


2 


G 
* = 0.005 + 0.011 = 0.016 
Vº 
v % 
= 0.126 2 — = 0.252 


V = (500 000 x 1.52) + 25.2º/ 
XxX 
ou V = (760 + 192) x 10º mº 


A maior participação do erro relativo no 
volume é dada pelo erro das espessuras minera- 
lizadas. Como poderemos calcular esse erro de 
um modo mais correcto? Veremos na segunda 
parte deste trabalho que o carácter estrutural 
(que é real) de uma mineralização pode ser 
introduzido pelo conhecimento do semivario- 
grama, função que permite calcular a variância 
de extensão da altura de cada poço ao seu 
rectângulo de influência de área ab. 


2. VARIÂNCIA DE EXTENSÃO E DE ESTI- 
MAÇÃO 


Vamos considerar a noção da variância de 
extensão como sendo a variância do erro que se 
comete atribuindo o valor de um dado parâme- 
tro mineiro (por exemplo, um teor) dum deter- 
minado domínio a outro domínio onde ele é, por 
hipótese, desconhecido. 
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Esse parâmetro será tratado como uma Variúvel 
Regionalizada, conceito introduzido por Matheron 
para designar aquelas variáveis cuja evolução no 
espaço, embora não obedecendo a nehuma lei 
rígida (conhecida), não pode ser considerada 
como puramente aleatória. É este carácter de 
regionalização que permite sujeitar a distribuição 
da variável a determinadas hipóteses, como vamos 
ver: Seja x a coordenada de espaço e Y (x) a 
variável em estudo. O facto de não podermos, 
sobre esta variável, aplicar a inferência estatis- 
tica (não podemos repetir provas no mesmo 
ponto do espaço, porque o teor num ponto é 
único) exige que consideremos o conjunto dos 
teores conhecidos como uma (e a única) realização 
de uma determinada Função Aleatória cuja lei 
de distribuição da probabilidade é evidente- 
mente impossível de determinar a partir de uma 
única realização. A hipótese que nos permite 
ficar a conhecer mais qualquer coisa sobre a lei 
de distribuição dos Y (x) é a chamada hipótese 
INTRÍNSECA que consiste em admitir a estacio- 
naridade dos quadrados dos acréscimos de Y (x) 
e que nos leva pois a definir uma função (x), 
chamada semi-variograma ou função intrínsica, que 
só depende do vector h que traduz o acréscimo: 


7 (h) = : E [Nov 6a | 


Desenvolvendo o quadrado e restrigindo um 
pouco mais a hipótese (introduzindo a estaciona- 
ridade de segunda ordem, que consiste em admitir 
que a variância e a covariância são independentes 
do ponto de apoio), temos que o primeiro e o ter- 
ceiro termos do desenvolvimento são iguais e a 
sua soma corresponde à variância à priori da 
função aleatória cuja realização é o conjunto dos 
valores conhecidos para Y (x). 

Nesta hipótese: 


v 4h) =K (0) — K (h) 


onde K (O) é a variância à priorie K (h é a 
covariância. 

Nas funções aleatórias usuais, K (h)tende para 
O quando h —» > (é intuitivo que uma medida 
da interdependência das duas variáveis se deve 
anular para uma distância infinitamente grande 
entre amostras) e então 7 (+) ==K (O). Assim 


(e sempre que houver uma variância à priori 
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finita para todos os h, por maiores que sejam, 
como é o caso na maioria dos jazigos de dimen- 
sões bem definidas, de fronteiras nítidas, espe- 
cialmente nos sedimentares (*), o semi-vario- 
grama será uma função crescente de h, até 
atingir um patamar cuja ordenada é K (O) que 
se inicia a uma abcissa b, dita amplitude, cujo 
significado físico é claro: é a distância a partir 
da qual passa a haver independência entre 
amostras. 

y (h) terá um andamento como o que se repre- 
senta na (Fig. 7), será uma função par e não 
negativa partindo da origem: 


Y th) 


K to) 


Fig. 7 


Evidentemente que haverá um semi-variograma 
para cada direcção do vector h (o que é um útil 
instrumento de análise e interpretação das ani- 
sotropias da regionalização). Como se constrói, 
na prática, um semi-variograma? Dispondo de 
uma carta do jazigo em que a variável regiona- 
lizada em estudo está referida às sondagens res- 
pectivas, escolhe-se a direcção « segundo a qual 
queremos calcular o semi-variograma e, segundo 
essa linha, vamos fazendo as diferenças X, —X, 
X.— Xo,... X,,+1— X, e dividimos os quadrados 
dessas diferenças pelo seu número n, obtendo 
assim o ponto ; (1). Fazendo agora as diferenças 
de dois em dois pontos, de três em três, etc, 
quadrando e achando a média dos quadrados, 
obtemos do mesmo modo y (2), y (3), etc. 
Como regra geral, pode dizer-se que o variograma 
não deve ter mais de 1/4 do número de pontos 
conhecidos para o jazigo. Assim obtemos um 


(*) Uma excepção a esta hipótese da estacionaridade 
de segunda ordem (mas que não tira generalidade a y (h) 
porque este é só definido com base na hipótese intrínseca), 
são os jazigos do Rand da África do Sul, em que Y th) 
aumenta indefinidamente com k (satisfazem ao chamado 
esquema de De Wijs—y (h) = 34 log h). 
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variograma prático ao qual assimilamos um 
«esquema» teórico (esquema linear, esférico, de 
De Wijs) e passamos a trabalhar sobre esse 
«esquema» que nos dá uma expressão analítica 
para y (h). Conhecido o semi-variograma, podemos 
iniciar sobre ele o cálculo das diversas variâncias 
de extensão e estimação. 

Designa-se por Variância de Extensão a variância 
do erro que se comete tomando o teor do volume 
v” (desconhecido) pelo teor do volume v (que se 
conhece). Designando por Z (v) e Z (v”) os teores 
médios nesses dois volumes, as suas expressões, 
admitindo uma distribuição uniforme da variável 
regionalizada são: 


Z (v) = [Y 60) dx 
V 
v 


zw) = [Y6) dx 
V 


Fig. 8 
A variância de extensão é 


2 A E aid 2 
cp=Eiz(v)—-z(v)p =0,—-20, +04 
onde o: designa a variância da variável regio- 


nalizada Z (v) em », qs , representa a mesma 
variância em v', e 7,,: , à covariância, 

Substituindo os teores médios pelas suas 
expressões, 


=5E).] Y (x) Y (y) dx dy 
V vd v 


Í 1 
= 2E[ [6 Y() dxady 
v 4 “ v' 


Ve 


1 
cp =2E [ | Y (y) Y (x) dx dy 
VV vo v* 
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Amitindo que se pode permutar o integral 
com a esperança (o que é válido porque a espe- 
rança corresponde a uma integração), 


= ( [ RE lic=s) de dy 
:, V wu Vv 


E = 31.) K (x—y) dx dy 
Vo dvd 


cg “fl K (x—y) dx dy 


co DRE 


onde K (x—y) é a covariância dos valores 
tomados pela variável regionalizada Y (x) em 
dois pontos quaisquer x e y (que podem varrer 
quaisquer dos dois domínios v e ví (ou os dois 
simultâneamente). 


Substituindo K (x-y) por K (0) —y (x - y) 
(o que corresponde, como vimos, a admitir esta- 
cionaridade de segunda ordem), obtém-se : 


o 
7 (x-y) dx dy — 
VV v' 


A [rtendranl a). | rtey) day 


que é a expressão geral da variância de extensão. 
Passemos agora ao caso descontínuo — o 


[4 


valor v é reconhecido por N sondagens pon- 


E 
tuais. Formemos um estimador Z (v), fazendo a 
média aritmética dos diferentes teores em v. 


é o verdadeiro teor 
do volume v 


Z (v)= poa À Y (y) dx 


* 1 da 
Z (v) = N » Y (x) é o estimador 


Designa-se por Variância de Estimação a variân- 
* 

cia do erro Z (v) — Z (v) e calcula-se como no 
caso anterior, obtendo-se: 


9 2 
= df te) da 


o A 
sh 7 | 7 (x- y) dx Hgroeçi 27 (4-9) 
v/v Nº 


q 


Como vemos, estas expressões dependem só do 
semi-variograma. Conhecido este, poderemos 
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calcular qualquer variância e portanto conhecer 
o erro que se comete estimando qualquer pará- 
metro através de um sistema de sondagens. 
Vejamos o caso prático mais interessante que 
é o caso do reconhecimento através de uma malha 
regular. Vamos ver como a fórmula geral atrás 
escrita se simplifica e se pode aplicar. 
Consideremos uma malha regular de sondagens 
de passo a: No centro de cada intervalo está 
implantada uma sondagem, como se representa 
na Fig. 9 (trata-se do chamado dispositivo centrado 
a uma dimensão). Vamos ver o erro que se 
comete estimando o teor médio do segmento L 
pela média aritmética dos teores das sondagens. 
Admitindo a independência dos erros, podemos 
supor que a variância de estimação que preten- 
demos é o quociente por N (número de sondagens) 
da variância de extensão elementar da sondagem 
implantada em x, ao seu segmento de influência a. 


Fazendo na fórmula geral da variância de 
estimação N==1 (uma sondagem) ei=1 
(o somatório tem um só termo), obtém-se notando 
quev=aeh=x—-y: 


= 
a 


* 


IE (x —y) dx — E ] | f (x —y) dx dy 


a/2 aa 
É caos mao E À E (ra 
=x = Jr ah cj y) dx dy 


Introduzindo as funções auxiliares 


h 


x (h) = à f bo) dx 


h h 


[fre dx dy 


oo 


1 


E (h) = 
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=. = —uue o e 


que estão tabuladas em (2), obtemos a expressão 
da variância de extensão para o dispositivo 
centrado a uma dimensão: 


o: =2%(a/)-—F(a) 


” 
Fig. 10 


No caso da Fig. 10, temos o chamado disposi- 
tivo fechado em que as sondagens estão implan- 
tadas nos nós da malha regular de passo a. 
Neste caso vem para a variância elementar: 


a 
9 


af 1 
o 7 (x— y) DP Pagar 
o 


1 - 

E | fro-y as dy 

fo 
a a 


ç? =2)1(9)-— 7 (a) — E (a) 


Procuremos aplicar as noções apresentadas ao 
cálculo da variância de extensão do erro que se 
comete estendendo o valor da espessura minera- 
lizada de uma sondagem à sua zona de influência. 

Para tal, já vimos que temos de adaptar um 
«esquema» teórico ao variograma experimental. 
Como o jazigo apresentado no ponto 1 é dema- 
siadamente pequeno para permitir uma adaptação 
correcta, trabalharemos sobre outro jazigo da 
mesma formação e suporemos que os resultados 
são extrapoláveis. O jazigo apresenta uma forma 
aproximadamente elíptica com duas direcções 
principais (correspondentes aos eixos principais 
da elipse) que designaremos por D; e Ds. 

Segundo D:, os valores das espessuras mine- 
ralizadas apresentam-se no quadro a seguir: 


i 
1 
2 0,70 
3 
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+ 0,00 
5 1,10 
6 1,70 
7 0,90 
8 3,40 
9 1,70 
10 3,30 
11 6,00 
12 0,55 
13 0,50 
14 1,00 
15 0,80 
16 3,50 
17 3,20 
18 1,50 
19 1,40 
20 3,80 
21 1,40 


O processo de cálculo já apresentado permitiu 
obter os valores seguintes para o variograma: 


27(1) = 3,5 
2/!2)==5,7 
27(3)=3,8 
27(4) =5,0 


Segundo a direcção perpendicular D,, os valo- 
res das espessuras mineralizadas são apresen- 
tadas no quadro seguinte : 


i se! 

1,30 
5,10 
2,95 
3,80 
1,30 
6,00 
4,50 
2,30 
4,80 
4,30 
1,50 
1,90 
1,30 


oO orn ãMa 


jd pd pá pa 
QnNm=Oo 


A existência de uma trincheira com amostras 
de 50 em 50m permitiu obter um ponto suple- 
mentar —y (1/2). 
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Os valores do variograma experimental segundo 
D, apresentam-se a seguir: 


2y(1/2)=1,25 


27(1) = 6,0 
27(2) =5,0 
27(3) =4,5 
27(4) =4,9 


Os valores segundo D; e Dy estão representa- 
dos no gráfico da Fig. 11 pelas linhas —.—.— 
e ——— , respectivamente. O cálculo da variân- 
cia à priori da distribuição dos X, segundo as 
duas linhas conduziu ao valor C=7(50) = 
=K (0) = 2,5. 

Uma interpolação gráfica aproximada realizada 
sobre as duas -linhas permitiu concluir que a 
relação entre as derivadas na origem segundo D, 
e segundo Dera 1,5: 1,0 que revela uma ani- 
sotropia geométrica que será destruída por uma 
afinidade de razão 1,5 aplicada sobre D, (a 
direcção de menor comprimento). Tal afinidade 
conduz ao traçado da nova linha corres- 
pondente a D,. Nesta altura, o jazigo é geome- 
tricamente isótropo. Se admitirmos que a tan- 
gente na origem que interpolará as duas linhas 
iscaais é a recta apresentada a cheio, 
esta intersepta a recta 2C == 5 num ponto (- 5/3) 
que, em «esquema» esférico de transição (a hipoó- 
tese de «esquema» conhecido mais adaptável aos 
resultados experimentais) corresponde a um ponto 
aí tal que a—3/2 à”. E assim fica definida a 
«amplitude» a. 

Se admitirmos então uma amplitude de 250 m 
e como a malha é de 100 m, podemos entrar no 


2Y (h) 


erercnencantes D2 depois da 
atinidade 


de módulo 1,5 
Fig. II 
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ábaco da Fig. 12 e encontrar na curva superior 
(para |/a = 0.4) uma percentagem da variância 
à priori de 15º/y. Concluímos portanto que, pelo 
facto de termos admitido a estrutura esférica para 
o jazigo (pelo facto de termos admitido aquele 
variograma), a variância à priori da distribuição 
das espessuras baixou a 0,15 X< 2,5. A variância 
do erro que se comete ao estender cada ponto à 
sua zona de influência é consideravelmente menor 
do que se tivessemos admitido (como fizemos no 
final da primeira parte) uma distribuição aleatória 
das espessuras mineralizadas. Se encontrarmos o 
variograma correspondente a uma determinada 
mineralização fica pois resolvido o problema 
deixado em aberto no fim do ponto 1 e os erros 
cometidos na estimação da altura e da espessura 
média são já da mesma ordem de grandeza. 


3. KRIGEAGEM 


Até agora atribuimos o teor da amostra central 
ao seu rectângulo de influência e, calculámes a 
variância do erro cometido (que é já inferior ao 
calculado na estatística clássica). O teor dessa 
amostra é um ESTIMADOR do teor do bloco. O 
problema que se põe naturalmente nesta altura 
é procurar o MELHOR estimador para esse teor, 
ou seja o estimador que minimize o erro obtido. 
Este é o princípio da krigeagem. Designa-se por 
krigeagem pontual o cálculo de teores em pontos 
desconhecidos, conhecendo o semi-variograma e 
teores em pontos próximos. 

Formulemos o problema: Um painel de teor 
Z desconhecido é amostrado por N sondagens 
de teores X,, X,... Xy. Em vez de dizermos 
que o teor do painel é a média aritmética dos 
teores amostrados e que a variância de extensão 


é a obtida pelas expressões do ponto 2 (onde Z 


é a média dos teores), vamos calcular os a, tais 


que o estimador Z === a, X, apresente variância 
Í 
mínima. 

A primeira condição é que o erro de Z — Z tenha 
esperança nula, o que dá = a,==1 (a soma dos 
factores de «peso» tem de ser igual à unidade). 

A condição de minimizar a variância do 
estimador formula-se derivando em ordem aos 

a 2 > 
a, a expressão de c (Z—Z). 
Essa expressão escreve-se: 


T(Z-D=0,-22a0,;+2aa;o [4] 
13 
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onde de N+1 equações a N-[-1 incógnitas (os a, e 
o À de Lagrange) que resolve o problema. 
ç; é a variância do teor real Z Como não interessa determinar ), vamos 
C,xi é à covariância da amostra X, com o teor procurar eliminá-lo. Toma-se uma sondagem 
do painel privilegiada (por exemplo no centro do painel) 
N—1 
%j é a covariância de duas amostras X, e X,. ax e faz-se ay = 1 — s á 
Tendo em conta a equação de ligação Z a, =1 Explicitando esse termo no somatório anterior, 
e aplicando a teoria dos multiplicadores de vem 
Lagrange, obtém-se o sistema N—l 
b a, 2, + dn Tin asi 2 ) 
e à, O = Co + À N—1 
2a =1 2 Pa PR a 
j 
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N—i 


eliminando 2, obtém-se N-1 equações > aR,—N, 
j= 


com (5) 


Destas N —1 equações, obtém-se os ã,, às «+ 
«ax. |, Sendo ay dado por diferença. 

Substituindo estes valores na expressão da 
variância já escrita (4) vem a chamada variância 
de Krigeagem: 


Esta é a variância minima que se pode obter 
(com a colecção de valores X,) para o estimador 


Z e tem uma interpretação curiosa: os três pri- 
meiros termos representam a variância de exten- 
são do teor de sondagem privilegiada X,, a todo 
o painel e, a essa variância, subtrai-se o soma- 
tório dos a, N, (termos sempre positivos) que 
vão baixar a variância tanto mais quanto maior 
for o número de pontos X;, de teor conhecido. 
As covariâncias dadas em (5) só dependem do 
variograma. 

A krigeagem pontual formula-se do seguinte 
modo : 

Conhecendo-se o teor em N pontos de coor- 
denadas x,, pretende-se conhecer o teor Z no 


EM 
ponto z. O estimador é da forma Z = 2a x 
e a resolução de um sistema semelhante ao ante- 


1 


rior conduz à mesma solução 


mas com 


R, =y(x,—x,)—7y(x,—x,)—7y(x, —x;) (6) 


N, =7(2 =x) -vbg=2) 764 =*) (7) 


Neste caso, os termos da matriz do sistema 
são lidos directamente sobre o semi-variograma, 
sem ser necessário calcular (a partir dele) cova- 
riâncias. 

Concretizemos a aplicação da Krigeagem pon- 
tual num exemplo concreto. Consideremos uma 
linha de sondagens num jazigo de ferro. 
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26.61 
29.80 
47.19 
52.53 
Xj > 57.84 
59.16 
61.63 
X;— 54.68 
X,— 57.18 
57.48 


60.46 
57.67 
48.92 
52.62 
57.02 
36.15 
52.66 
48.51 
91.27 
58.54 


Os teores de Fe do «tal-qual» obtidos em son- 
dagens ao longo dessa linha através de uma 
malha regular de sondagens de passo de 100 m, 
estão representados no quadro ao lado. Suponha- 
mos que não conhecíamos o teor, por exemplo 
no ponto 7. Vamos calculá-lo por krigeagem 
(ponderando adequadamente os teores nos outros 
pontos, por hipótese conhecidos), determinando 
depois a variância de Krigeagem (a menor variância 
possível para esta colecção de valores) e verifi- 
cando que o teor no ponto «desconhecido» está 
contido no intervalo de confiança dado por essa 
variância. 

O primeiro passo a dar será encontrar um 
«esquema» e adaptá-lo ao variograma experi- 
mental. Os valores 2 y (h) estão representados 
na Fig. 13 em função de . | 

O teor médio é LX,/19=50.80e 5X; /19 = 
=— 2680, o que dá para a variância, 100 pontos 
quadrados. O patamar no variograma (como 
estamos a usar 2% (h)) será 200 e a análise da 
distribuição dos pontos experimentais sugere a 
aplicação do «esquema» linear. Este variograma 
tem a particularidade de apresentar o chamado 
«efeito de pepita» (ordenada na origem da recta), 
designado por 2C, e que corresponde ao facto 
de não termos valores do variograma para |: 
inferior a 100 m, 
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220 — 


200 -—+ | | | 


180 — 


Então teremos: 


2 €, (efeito de pepita) = 33 
2 C (patamar) == 167 
a (amplitude) = 450 


A equação da recta (reduzida ao esquema linear 
teórico) será: 


2c 167 


2y(h)= — = —— para h<450m 
a 450 


2y(h)== 2c = 167 para h > 450m 

As equações da krigeagem serão representadas 
por uma matriz 18 >< 18 que nos dão os a,. Como 
esta matriz contém muitos termos iguais (surgi- 
mento do «efeito de ecran» — os teores mais pró- 
ximos do ponto a Krigear mascaram a influência 
dos teores mais longínquos) tentámos a resolu- 
ção das equações da Krigeagem usando apenas 
os pontos x,, X,, Xy e x, que rodeiam o ponto 
a krigear (quadro anterior) e considerámos o ponto 
x, como um ponto privilegiado situado no centro 
do jazigo e com o teor médio (50.80) de todos 
os pontos. Então, designando por 2C/a o coefi- 
ciente angular da recta podemos escrever, de 
acordo com a expressão geral (6), e atendendo 
a que a matriz é simétrica: 
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2C 7167 


q 0 —— sm e ——— e—t 
——— — “ + + 4 +. + 


h.1072 (m) 


Ru=—27(55)=—2% 


Ro=7(D—y(4.5)—y(5.5)=—2e+ — 


Ro =—27(4.5)=—2c 
c 5€c 
R5=-c+t— Ryg=-c+— 
a Za aê 2a 
c 3c 
R;=-c—— Ry=-c+— 
25 
2 =” 2a 
5c 3c 
Ra = — — Ray = — — 
a a 
Sc 
Rag e 
Do mesmo modo, as equações [7] dão: 
N=-e-2€ 
Za 
5 
N=-c— > = 138.2 
a 
E 
N=—-— =4.65 
a 
3c 
N, = — 
a 


a matriz reduz-se a: 


98427] [a |=767 
8953 a, | 
24553 | | a, 5 | 
2783] | a | 


a 


415 


que dá a=a=0 


3, = à = 0.5 


Verificamos que neste caso o efeito de ecran 
é total e que os dois pontos mais próximos do 
ponto a krigear são, por si sós, responsáveis pelo 
teor no ponto desconhecido. 

Então esse teor (obtido por krigeagem) será 
56.97 e o verdadeiro valor é 61.83. Vejamos se 
está compreendido no intervalo de confiança da 
variância de krigeagem. A expressão desta 
variância é (no caso da krigeagem pontual) : 


=27(Z-x,)+ZanN, 


ou seja 27 (3,5) — 0,5 138,2 — 0,5 X< 46,5 = 
= 13,4, o que dá um intervalo de confiança do 
dobro do desvio padrão de 7,32, que compreende o 
verdadeiro valor do teor no ponto «desconhecido». 
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J. J. D. Domingos C. D. U. 532.5-4 + 532.54 


Recirculating flows 


Técnica No. 407 — XLVI — 4, 1971, pág. 347-356. 


Discutem-se os escoamentos bidimensionais em que as 
aproximações de camada limite não são válidas, e os 
aspectos matemáticos e físicos do problema, As dificul- 
dades intrinsecas do tratamento rumérico das equações 
são exemplificadas, mostrando-se que o aumento do nú- 
mero de Reynolds conduz a exigências contraditórias 
no que se refere a precisão e estabilidade, 
Apresentam-se modelos de turbulência e a sua incorpo- 
ração num método de cálculo aplicável a este tipo de 
escoamentos tal como desenvolvido por D. B. Spalding 
et al. 


S. D. Antunes C.D.U. 551.508.5 : 53.089 


An introduction to hot wire anemometry 
Técnica No. 407 - XLVI— 4. 1974, pág. 357-366. 


Justifica-se o uso de fios quentes em anemometria, em 
particular descreve-se o método de calibração destes, 
e interpreta-se a resposta do anemómetro de fio quente 
em escoamentos de baixa turbulência. 


Duarte - Ramos C. D. U. 519.738 
Sobre o critério Liapunov 
na teoria da estabilidade de sistemas 


Técnica No. 407 — XLVI — 4. 4974, pág. 367-379. 


Definem-se as noções de estabilidade e efectua-se a sua 
interpretação geométrica. Exemplifica-se o método directo 
e o critério Liapunov num sistema mecânico. Compa- 
ram-se os diversos métodos de formação de funções Lija- 
punov em sistemas não lineares, 
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Influência das tensões negativas nas características 
estruturais dos solos compactados 


Técnica No. 407 — XLVI — 3. 1974, pág. 381-350. 


O presente trabalho refere-se ao estudo de solos argilosos 
compactados, e o seu objectivo imediato é a aquisição de 
elementos que possam contribuir para o conhecimento 
da reologia desses solos em termos de tensões efectivas. 
Foi dada especial atenção à medição das tensões neutras 
instaladas na fase liquida do solo parcialmente saturado, 
bem como às técnicas que permitem conhecer a evolução 
dessas tensões durante os ensaios de determinação de 
caracteristicas mecânicas. 

Avaliou-se ainda experimentalmente, a influência do 
estado de tensão endógeno, associado à situação de satu- 
ração parcial do solo, no comportamento desse mesmo 
solo após a saturação, 

Os resultados indicam que o grau de saturação de fabrico 
determina todo o posterior comportamento reológico dos 
provetes. 
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Aválise de estruturas redutíveis 
a osciladores não lineares com histerese 


Técnica No. 407 — XLVI — 4. 1971, pág. 391-402. 


Nos estudos sobre Engenharia Sismica que vêm sendo 
conduzidos no LNEC desde há cerca de 15 anos, tem 
sido necessário aplicar técnicas de anúlise numérica 
a osciladores de um e vários graus de liberdade, linea- 
res e não lineares, os quais constituem os modelos mate- 
máticos de diversos tipos de estruturas. ; 
Descrevem-se neste trabalho as características dos osci- 
ladores que têm sido analisados, relacionando-as com 
os tipos de estruturas que representam. Referem-se os 
principais aspectos da formulação das equações diteren- 
ciais que descrevem o comportamento dos sistemas e das 
técnicas utilizadas na respectiva integração «passo a 
passo», Apresentam-se aplicações ao estudo do compor- 
tamento de edifícios de betão armado actuados por Sis- 
mos de forte intensidade, 
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Aplicações da geoestatística 
Técnica No. 407 — XLVI— 4. 1971, pág. 4032-416. 


Em três problemas práticos, pretende-se provar que os 
métodos da geoestatística, na medida em que conferem 
ao depósito mineral uma certa «estrutura» que lhe é dada 
pelo modo particular como ocorreu a mineralização, con- 
duzem a menores erros do que um tratamento pela esta- 
tistica clássica. 
Esses problemas são : 

|. Definição da superficie e volume mineralizado. Erro 
deste cálculo, 

2. Variância de extensão de uma amostra à sua zona 
de influência. 

à. A Krideadem — ou seja a procura do melhor estimador 
para um certo parâmetro (aquele que minimize a variân- 
cias definida em 2,). 
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Mananua DAS NEvES 


The influence of negativo neutral stresses on the 
structural characteristics of compacted soils 


Técnica No. 407 — XLVI — 4. 4974, pag. 3814-590, 


The present paper studies compeacted clay soils in order 
to obtain experimental data which may help to investigate 
their rheolody im terms of effective stresses, 

Two points specially considered were the measurement of 
the neutral stresses in the liquid phase of a partly satu- 
rated soil and the procedures for determining the evolution 
of these stresses during tests ofits mechanical properties. 
The infinence on the behavionr of the soil after satura- 
tion, of thz endogenijc state of stress due to the initial 
partial saturation, was siso experimentallyv studied. 
Accord ng to the results, the tabric percent saturation 
entirely determines the subsequent rheologic behaviour 
of the specimens, 
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Artur Ravana 
MADALENA (QuIRIRO 
Mania Tenesa Peneira DA SILVA 


Dynamic analysis of non linear hysteretic systems 
Técnica No. 407 — XLVI — 4.1971, pag. 3914-402. 


In the studies on earthquake engincering that have been 
carried out at the Laboratório Nacional de Engenharia 
Civil for about 15 vears, different techniques have been 
applied for the analysis of single and multidegree of 
freedom oscillators, linear and non-linear. 

The present paper describes the main characteristics vf 
the systems that have been considered, relatiny them to 
different types of structural behaviour. The formulation of 
the differential equations of motion as well as the step- 
-by-step integration techniques are briefly reported. 
Examples of dynamic analysis of reinforced concrete 
structures acted by strong errthquakes are presented. 


Garcia PenremA 
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Geostatistic Applications 


Técnica No. 407 — XLVI — 4.14971, pag. 403-416. 


In three pratical problems, it is intented to prove how 
geostatistic methods,as theytakeinaccountthe «structure» 
given to the ore deposit by its particular mineralization, 
lead to bether significance levels (to smaller errors) than 
classical statístics. Those problems sre: 

1) Definition of the deposit surface and volume. Error 
of th's calculation. 

q, Extension variance of a sample to its influence zone. 

») The Krigeadge — the search of the best estimation of 
a certain paramether minimizing the variance treated 
in 2.). 
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Recirculating flows 
Técnica No. 407 — XLVI — 4.4974, pag. 347-356. 


Two dimensional flows where the bcundary-laver appro- 
ximations are not ecceptable are introduced and their 
specific mathematrica! end physical problems discussed. 
Inerent difficulties to the numerical treatmentofthe partia] 
differential-equations are exemplified and it is shown hcw 
increasing the Revnolds number imposes contradictory 
requirements regurding stability and accuracy, 
Turbulence models are imtroduced together with a 
numerical solution procedure appropriate to these type 
of fiows, as developed by D. B. Spalding et al. 


S. D. Antunes U. D. C. 551.508.5: 53.089 


An introduction to hot wire anemometry 
Técnica No, 407 — XLVI — 4. 1971, pág. 357-366. 


This paper provides an introduction to the physics of hot 
wire anemometry. The use of wires as transducers for 
turbulence measurements is justified, and for flows with 
low turtulence the means of signal interpretation is 
discussed. 

Techniques of calibration are also described, In parti- 
cular, the constant temperature (constant resistance) and 
constant current hot wire anemometers are presented, 
and their uses are commented on. 
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H. Duanre Ramos U. D. CG. 519.738 


Liapunov criterion on stability of systems 


Técnica No. 407 — XLVI — 4. 1971, pág. 367-379. 


We define the conceptions of stability and show their 
geometric interpretation. The direct method and the 
Liapunov criterion are applied on a mechanical system. 
The methods of formation of Liapunov functions in non- 
“linear systems are compared. 
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